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ABSTRACT 


This  report  formulates  the  equations  of  motion  for  a  gun  tube  with  arbitary  space 
curvature  2md  cross  section  geometry  and  presents  results  showing  the  effects  of  curvature, 
torsion,  and  cross  section  shape  on  natural  frequencies  and  mode  shapes. 

Expressions  for  the  kinetic  and  potential  energy  functions  of  a  curved  and  twisted 
tube  are  formulated  in  terms  of  centerline  displacements  2md  cross  section  rotations.  These 
expressions  are  first  used  to  derive  the  equation  of  motion  and  natural  boundary  conditions 
using  Hamilton's  Principle.  The  equations  of  motion  are  then  specialized  for  various 
variations  of  centerline  curvature,  torsion  and  cross  section  shape  along  the  tube  axis. 

Next,  the  kinetic  and  potential  energy  functions  are  used  to  formulate  a 
Rayleigh-Ritz  solution  for  the  mode  shapes  and  natural  frequencies  of  helical  tubes  with 
fixed-fixed  ends.  The  effects  of  space  curvature  and  torsion  and  of  cross  section 
out-of-roundness  on  frequencies  and  mode  shapes  are  examined. 

The  natural  frequencies  and  mode  shapes  of  tubes  are  found  to  be  significantly 
affected  by  even  a  slight  change  in  values  of  the  curvature  and  torsion  parameters.  The 
fundamental  frequencies  and  the  degree  of  coupling  among  the  displacement  variables  in 
natural  modes  always  increased  when  either  curvature  parameter  was  increased. 
Out-of-ro’:.  Jness  edso  introduces  further  complexity  into  the  mode  shapes. 


X 


CHAPTER  1 


INTRODUCTION 

The  work  presented  in  this  report  was  sponsored  by  the  United  States  Army  for 
the  purpose  of  examining  the  effects  of  space  curvature  and  cross  section  out-of-roundness 
on  the  vibration  mode  shapes  and  natural  frequencies  of  gun  tubes. 

In  the  context  of  structural  mechanics,  the  problem  is  that  of  vibrations  of  curved 
and  twisted  rods  of  arbitrary  cross  section.  The  problem  of  vibration  of  rods  curved  in  a 
plane,  especially  circular  arcs,  has  been  the  subject  of  a  large  number  of  investigations.  It  is 
found  that  this  type  of  curvature  results  in  mode  shapes  in  which  lateral  or  out-of-plane 
motion  is  coupled  with  twisting  and  in-plane  motion  in  which  in-plane  axial  and  transverse 
displacements  are  coupled.  The  effect  of  three  dimensional  curvature  on  rod  vibrations 
does  not  appear  to  have  been  extensively  studied.  Indeed,  no  paper  on  this  subject  could 
be  found  in  a  survey  of  the  vibration  literature  of  the  last  twenty  years. 

Although  the  vibration  theory  of  curved  rods  has  great  importance  in  many 
engineering  applications  and  has  attracted  the  interest  of  many  investigators,  a  completely 
general  formulation  which  might  illuminate  unnoticed  dynamic  effects  and  interactions 
caused  by  such  phenomena  as  initial  space  curvature,  variable  shape  and  area  of  cross 
section,  and  initial  twist  of  cross  sectional  principal  coordinates  along  the  rod  axis  is 
difficult  to  find  in  the  literature. 
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One  of  the  seminal  treatments  of  the  derivation  of  equations  on  the  motion  of  rods 
with  initial  curvature  appears  in  the  text  on  elasticity  by  Love  10  .  Although  restrictive 
forms  of  Love's  equations  are  used  directly  or  are  re-derived  by  subsequent  investigators  for 
small  amplitude  vibration  of  rods  curved  in  circular  arcs,  there  are  deficiencies  in  Love’s 
work  which  render  his  strain-displacement  equation  unsuitable  as  the  starting  point  for  a 
general  examination  of  the  motion  of  rods  with  space  curvature.  In  Love’s  work,  the  final 
state  of  deformation  is  assumed  to  be  such  that  cross  sections  remain  plane  euid  normal  to 
the  center  line  of  the  deformed  rod.  By  these  assumptions,  not  only  is  transverse  shear 
deformation  excluded  but  also  kinematic  inconsistencies  are  introduced. 

Love's  displacement  was  recently  re-examined  by  Kingsbury  8. 9  who 
formulated  governing  equations  for  a  curved  and  twisted  rod.  Kingsbury  employed  a 
completely  general  set  of  small  strain  but  large  displacement  and  rotation 
strain-displacement  equations  with  no  a  priori  assumptions  regarding  the  structural  action 
or  mode  of  deformation  of  the  rod.  Thus,  the  functional  form  of  the  displacement 
components  was  unspecified.  A  technical  theory  of  rods  with  space  curvature  was  then 
developed  based  on  the  assumption  that  cross  sections  remain  pleme  and  undeformed.  A 
system  of  four  coupled  equilibrium  equations  in  four  displacement  and  rotation  variables 
resulted. 


Kingsbury's  technical  theory  of  curved  rods  provided  the  starting  point  for  this 
study.  When  a  normal  mode  solution  for  the  free  vibration  of  helical  rods  was  attempted, 
however,  it  was  found  that  natural  frequencies  could  not  always  be  obtained  from  the 
resulting  characteristic  equations.  It  was  felt  that  this  problem  could  be  the  result  of  terms 
in  the  governing  equations  which  caused  the  stiffness  matrix  of  the  characteristic  equation 
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to  be  non-symmetric.  A  second  problem  encountered  was  that  of  uncertainty  concerning 
the  proper  form  of  the  force  boundary  conditions  for  the  rod.  It  was  finally  realized  that 
approximate  solutions  for  mode  shapes  and  natural  frequencies  for  a  variety  of  boundary 
conditions  could  most  easily  be  obtained  by  using  the  Rayleigh-Ritz  method  which  requires 
the  use  of  potential  and  kinetic  energy  functions.  For  these  reasons,  it  was  decided  to 
undertake  a  new  formulation  of  the  theory  using  Hamilton’s  Principle  or  the  Principle  of 
Virtual  Work  for  dynamical  systems. 

The  derivation  of  the  governing  differentizd  equations  employs  the 
strain-displacement  relations  introduced  by  Kingsbury  '8,  9  .  These  strain-displacement 
equations  are  based  on  the  assumptions  that  displacements  are  small  and  the  radius  of  the 
rod  cross  section  is  small  in  comparison  to  the  radius  of  curvature.  In  addition,  it  is 
assumed  that  the  cross  sections  remain  normal  to  the  centroidal  curve  locally  (at  centroidal 
axes),  or  in  the  sense  of  vanishing  of  the  average  transverse  shear  strain.  The  formulation  of 
strain  energy  retains  the  effects  of  extension  of  the  centroidal  curve. 

In  Chapter  3.  Hamilton's  principle  is  used  to  formulate  the  governing  equations 
and  natural  boundary  conditions  for  a  curved  rod.  A  system  of  four  partial  differential 
equations  in  four  displacement  variables  is  obtained. 

In  Chapter  4.  the  governing  differential  equations  (Euler-Lagrange  equations)  are 
expressed  in  terms  of  displacements,  cross  section  properties,  centroidal  axis,  torsion, 
curvature,  and  their  derivatives  for  specizd  cases. 


In  Chapter  5.  the  strain  energy  function  for  the  helical  rod  is  employed  with  the 
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Rayleigh-Ritz  method  to  formulate  approximate  solutions  for  mode  shapes  and  natural 
frequencies  for  helical  rods  with  both  ends  fixed  (clamped-clamped)  and  with  one  end  fixed 
and  the  other  end  free  (clamped-free).  Two  different  trizil  functions  (admissible  functions) 
are  used.  Use  of  both  types  of  functions  allows  comparison  of  frequencies  zmd  mode  shape 
results  obtained  from  those  functions. 

The  equations  of  motion  of  curved  rods  derived  in  this  report  are  compared  with 
those  derived  by  Kingsbury  8.  9  and  slight  differences  are  noted  between  the  two.  The 
governing  equations  of  a  ring  are  then  compared  with  those  presented  by  Love.  Chapter  6 
assesses  the  quantitative  effect  of  these  differences  by  employing  the  simplified  strziin 
energy  functions  with  a  corresponding  one  presented  by  Love  to  obtain  Rayleigh-Ritz 
solutions  for  natural  frequencies  of  the  out-of-plane  vibrations  of  a  clamped-clamped  ring. 

In  Chapter  7.  the  Rayleigh-Ritz  formulation  for  clamped-free  helical  rods  shown 
in  Chapter  5  is  used  to  examine  the  effects  of  curvature,  torsion,  and  out-of-roundness  on 
the  fundamental  natural  frequencies  and  mode  shapes  of  a  helical  rod. 

Finally.  Chapter  8  briefly  summarizes  the  conclusions  drawn  from  the  analyses 
performed  in  the  course  of  this  report. 


CHAPTER  2 


STRAIN  ENERGY  OF  A  GENERAL 
CURVED  AND  TWISTED  ROD 

In  this  chapter  the  strain  energy  function  for  a  rod  with  general  space  curvature 
and  cross  section  geometry  is  developed.  The  strain  energy  is  phrased  in  terms  of  cross 
section  displacement,  rotation  components,  and  the  curvature  and  torsion  of  the  rod 
centroidal  axis. 


2.1  Linearized  Strain-Displacement  Equations 


Linearized  strain-displacement  equations  for  a  general  curved  and  twisted  rod  are 
presented  in  a  paper  by  Kingsbury  8.  9 .  These  equations  were  based  on  the  assumptions 
employed  in  the  strength  of  material  formulations  for  the  bending  and  torsion  of  prismatic 
rods.  It  was  assumed  that  a  typical  cross  section  would  be  unchanged  in  size  and  shape 
after  the  structure  was  deformed,  and  each  cross  section  w’as  assumed  to  undergo  a  small 
rotation  around  each  of  the  three  coordinate  axes  while  its  centroid  was  displaced  along 
each  of  these  axes;  the  resulting  strain-displacement  equations  were  : 


d  w  da,.  da,  ,  .  . 

€„  =  _-  1-^-  y-^- K(u-  y<t>) 
d  z  a  z  d  z 


2c„  =  ^-Q-Av-xu;-y(KQ,-^)-z(Ad>-xQ) 
dz  d  z 


(2.1) 

(2.2) 

(2.3) 
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c 


24^  =  1-^+ A  U  +  Q,+ (2.4) 
d  z  a  z 

Equations  (2.1),  (2.2),  (2.3)  and  (2.4)  are  called  the  strain-displacement 
equations.  In  Elquations  (2.2),  (2.3)  and  (2.4),  u,  v  and  w  are  the  displacement  components 
at  a  point  on  the  centroidal  curve  (the  origin  of  the  local  x,  y  and  t  coordinates);  A  and  k 
are  the  initial  torsion  and  curvature,^  while  o^  o^  and  are  sntall  angles  of  rotation  of  the 
cross  sections  about  the  x,  y  and  z  axes  respectively  as  shown  in  Figure  2>1 


Figure  3>1:  Displacement  and  Rotation  Components 
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It  is  noted  that  as  in  the  technical  theory  of  curved  beams,  the  non-linear 
functions  and  the  additional  terms  representing  the  warping  of  cross  sections  of 
non-circular  rods  are  not  considered. 

U'hen  the  cross  sections  are  assumed  to  remain  normal  to  the  centroidal  curve 
locally  at  I  =  0  and  ^  =  0.  or  in  the  average  sense  by  integrating  Equations  (2.3)  and  (2.4) 
over  the  cross  section,  the  equations  expressing  the  vanishing  of  transverse  shear 
deformation  result  in  the  following  expressions  for  the  cross  section  rotations  : 


Q 


y 


d  u 
d  z 


X  u  —  K  w 


(2.5) 


a 


I 


d  V 

d  z 


A  u 


(2.6) 


2.2  Strain  Energy  of  a  General  Curved  and  Twisted  Rod  with  Doubly 
Symmetrical  Cross  Section 


The  strain-displacement  equations  introduced  in  the  last  section  will  be  used  for 
the  derivation  of  the  strain  energy  function.  Assuming  the  material  constituting  the  rod  to 
be  linear  elastic,  the  strain  energy  for  the  rod  can  be  expressed  by  : 


1  / 

2 


U=-  /  (^(7dV 


The  stress  and  strain  vectors  in  the  above  expression  are  : 


’’ly  ’ '  yj  •  . 


•  ~1iy  ’  lyi  ’  Ixi . 


(2.7) 


(2.8) 

(2.9) 


The  applicable  stress-strain  relations  are  : 
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'it  -  2  G€„ 


(2.10) 

(2.11) 

(2.12) 


where  E  is  the  modulus  of  elasticity  and  G  is  the  modulus  of  elasticity  in  shear. 


Upon  substituting  Equations  (2.1)  through  (2.6)  and  Equations  (2.8)  through 
(2.12)  into  Equation  (2.7).  the  following  expression  for  the  strain  energy  is  obtained  : 


U-  /  Piiy-P2ii-P3yj/-P4-P5i-P6VdV’ 

J  V 

where  the  quantities  P j.  Pj.  P3.  and  P^  are  presented  in  Appendix  A. 


(2.13) 


Since  the  centroidal  axis  is  the  origin  of  the  coordinate  system. 


xdA  =  0 


L 

I  ydA-O 

J  A 


(2.14) 


(2.15) 


Pj  and  Pg  terms  vanish.  Hence,  the  strain  energy  is  simplified  to  the  following  expression  : 


U=  f  Pi  I  -Pj/  -P3/  -P4A  dr 

y  *  ly  i  X3  J  yj-  * 


(2.16) 


where  I  is  the  length  of  the  rod.  According  to  Equation  (2.16),  the  strain  energy  function 
per  unit  length  of  rod  is  given  by  : 


S.E.F.  =  P,I  -P, ;  -Pol  -Pi  A 

*  ly  *  B  yy  * 


(2.17) 


where 


9 


4= 

f  yUA 

(2.18) 

J  A 

•^yy  ~ 

f  dA 

(2.19) 

J  A 

/  xydA 

(2.20) 

^A 

f  dA 

A 

(2.21) 

It  is  noted  that  the  area  moment  inertias,  area,  torsion,  curvature,  and  the  tensile 
and  shear  moduli  may  arbitrarily  vary  with  z.  the  coordinate  along  the  curve.  Although  the 
average  transverse  shear  strain  has  been  set  equal  to  zero  by  the  use  of  Equations  (2.5)  and 
(2.6).  extension  of  the  centroidal  curve  is  retained. 


CHAPTER  3 


derivation  of  the  governing  equations 

FOR  A  GENERAL  CURVED  AND  TWISTED  ROD 

In  this  chapter  the  differential  equations  governing  the  deformation  of  the  general 
curved  and  twisted  rod  are  developed  by  using  Hamilton's  principle.  For  a  continuous 
system,  the  use  of  Hamilton's  principle  provides  an  important  advantage;  namely,  the 
natural  boundary  conditions  are  found  systematically  in  the  process  of  obtaining  the 
equations  of  motion.  Knowledge  of  the  boundary  conditions  is  necessary  for  determining 
natural  frequencies  by  the  Rayleigh-Ritz  method,  which  will  be  discussed  in  a  later 
chapter. 

3.1  Hamilton’s  Principle  Applied  in  a  Curved  and  Twisted  Rod 

Hamilton's  principle  may  be  stated  as  : 

^6W  dt  =  0  (3.1) 

nc 

where 

T  =  total  kinetic  energy  of  the  system. 

U  =  potential  energy  of  the  sj'stem,  including  the  strain  energy  and  the 

potential  energy  of  conservative  external  forces. 

=  virtual  work  done  by  nonconservative  forces,  including  damping  forces 

flC 
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and  external  forces  not  accounted  for  in  U. 

b  =  variational  operator. 

The  potential  energy  and  kinetic  energy  functions  are  given  by  : 

(3.2) 

(3.3) 


and  p  —  mass  density. 


Jo 

T=  f\p  Ij0  dz 

Jq  2 

where 


According  to  the  Equation  (3.3).  the  kinetic  energy  function  is  defined  by  : 

K.E.F.=  -^p  A(ti'2-f-2-u-2)-(4-/y^)o-2:  (3.4) 

In  Equation  (3.2)  and  Equation  (3.3).  u.  r.  w.  o  and  their  derivatives  can  be  considered  to 
be  functions  of  the  independent  variables  zand  t. 

For  the  undamped  free  vibration  problem,  there  are  no  nonconservative  or 
external  forces  acting  on  the  system.  Combining  Equations  (3.2)  and  (3.3),  substituting 
into  Equation  (3.1),  then  integrating  by  parts  and  equating  the  resulteuit  coefficients  ofiu. 
Sv,  6w  and  So  yields  the  following  Euler-Lagrange  equations  as  well  as  the  natural  boundeu'y 
conditions  for  the  system. 

a*  dF  _  ^dF  _dF  _  d  dF 
a?a^7^  dzdu  d  v!  dtdu 


(3.5) 
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dF  d  dF  dF  d  dF 


d  ^  d  dzdv  dv  dtdv 
d  dF  dF  d  dF 


=  0 


d  zd  w  d  w  d  td  w 


=  0 


where 


^dF  dF_  d  dF 
d  zd  o'  do  dtdo' 


d  z 


_  in.  in.  I  .  I  . 

J-  I  U  .  U  U  .  U  .  V  .  t  .  f  .  V  .  w .  u-  .  w  .o  .o  .  u 


-p  A  I  -  v'^  ~  ~  ly^)o' 


(  Pi  ~  ^2  ^rr  '  ^3  ^yv  “  *^3  ^ 


and  a:  r  =  U  and  r  =  I 


Either  u  specified  or  -  - — j-j  —  0 

d  ti 


Either  u  j,«  specified  or - ^  - - j  =  U 

d  zd  u  d  u 


Either  v'  is  specified  or  ^  =  0 

d  V 

r*k  t  j  d  dF  dF 

Either  v  w  specified  or - ^  - - j  =  0 

d  zd  V  d  V 

d  F 

Either  w  is  specified  or - j  =  0 

d  w 

d  F 

Either  o  is  specified  or - ,  =  0 

d  o 


(3.6) 

(3.7) 

(3.8) 


(3.9) 


CHAPTER  4 


EXPRESSION  OF  EULER-LAGRANGE  EQUATIONS 
FOR  GENERAL  AND  SPECIAL  CASES 

In  this  chapter,  the  Euler-Lagreinge  equations  for  a  rod  with  space  curvature  are 
expressed  in  terms  of  the  rod  displacement  functions,  torsion,  curvature,  section  properties, 
and  their  derivatives.  The  general  equations  are  then  specialized  for  various  cases  of 
centerline  curvature. 

4.1  General  Form 

Upon  substituting  the  function  F.  Equation  (3.9).  into  Equations  (3.5).  (3.6).  (3.7) 
and  |3.  j).  the  following  equations  are  obtained  ; 

rr 

Ely^U- 

-(3EI^'-{2EI^'-EXI^)k)w" 

-{ZEl^  "-[\ElJ~  2EX1^)k  G1^^-(EI^"-EXJ^'-  AE)k]w' 

-{Eiy-{2ElJ-EXI^)K" 

-(-ZGI^^~EiJ'~EXI^')k'-{GXI^-GiJ)k^)w 

FT  "" 

-EI^V 

-(-ea/^-2£;^'-£a;Jv"' 
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-{~Gl^^-2EXlJ-2EX'ly^-ElJ'-EX^I^-EXlJ)v" 
-(-  2GI^k'-{GX1^-  GI^'~GXIJk^ 

-  EXlJ 4EX  'ij- 3£A ' ' I^-EX^lJ ~2EXX ' l^)v' 
-{2GX1^k'~{GX1J-Gx'i^-GX\)k^ 

-  EX  'ij 2EX '  'ij- EX ' ' 'l^~EXX 'iJ-EXX '  'l^)v 

r-T  "" 

-EL..U 

yy 

-2ElJu"' 

-i-GI^^-ElJ'-EXI^'-ZEX'l^-EX^IJu" 

-{-2GIy^K'-GljK^ 

- EXl^ '-4EX'lJ-ZEX ' 'l^- EX^lJ- 2EXX 'lju' 

-(-2GXI^k' GXlJ -  GX' I^-GX'^I^-AE)k'^~EX' 
-2EX  '  'l^'-EX ' '  EXX  'ij- EXX  " Iju 

-(-  G-  E)l^o" 

-({-G-2E)iy-{-GXl^-{-G-2E)lJ-{G-E)XIjK)o' 

-  (-  £ V  " -  GX/^-  2EI^'~EXI^)k  ' 

GXlJ-  GX' EI^"-GX^I^+EXiJ)k)^  =  -  Apt*" 

E>r  in 

EI^w 

-iZEI^'-{EXI^-2El^')K)w" 


(4.1) 
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-{E1^"'~(EXI^-2EiJ)k" 

-(-  ZG1^^-EXIJ-E1^")k'-{-  gxi^-  giJ)^)w 
-ElJ'" 

-2ElJv" 

- (-  G/^- 2-  £a2;^-  EXiJ- ZEX 'l^-ElJ'W ' 

-  (-  2GI^k  '- GIJk^~  EX^iJ-  2EXX  '/^ 
-EXlJ'-AEx'l^'-3Ex"l^)v' 

-(2GXI^k'-{GX^I^~GXI^'-GX'i^)k^~EXX'iJ 

-EXX'  'l^-  EX  2EX '  'j^'- EX '  "l^)v 

rr  "" 

-[EXI^-2EI^'-EXI^U"' 

-(-  GI^^-EXlJ -ElJ' ~EX^I^-2EXlJ^ZEX' IjJ' 
-(-  2G1^k  '-(-  GXl^-  GlJ- GXIJk^-EX^I^' 

-2EXX  'l^-EXlJ'-AEX  'lJ-ZEX"lJu 

-(-  2GXJ^k  GX%-  GXlJ-  GX  'iJk^-EXX 

^EXX '  'l^-EX  'lJ'-2EX '  'iJ-EX ' '  '/Jti 
-{-G-E}I^o" 
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-[[-  G-2E)I^' -[{-2G~  E)\I^-{-  G-2E)iJ)k)o' 

- (-  G-  E)X1^-  2EIJ)k  ' 

-(-  GAV^-I-  G-£:)A/^'-GA'/^-£;/„' V)o  =  -Apv"  (4.2) 

{-EI^^-AE)w" 

-(-2EI^k'-EIJk^-a'e)w' 

-{-EI^k"-EiJkk'-GI^*)w 

rr 

-(EM^-EjJ)kv" 

-{GI^^-(EXiJ-2EX'i^)k)v’ 

-{{Ex'iJ-E\"i^)k- GX1^^)v 

-i-ElJ-EXI^Wu" 

-{GI^^-{- EXiJ-  2Ex'i^-AE)k)u' 

~(AEk'  ~GXI^^-{-  ex'  I^'- ex"  I^-A'  E)k)u 
-{G-E)I^^o' 

-{EI^<'-{GXJ^-EI^')k^)<P  = -Apw"  (4.3) 

{-G-E)I^^w' 
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-((-  2G-  E)1^k  GI^')k'^)w 

-{-G-E)I^v" 

-(i(2G-E)XI^-GjJ)K-GI^')v 

-(GX1^'-(-GX^J^-GXIJ-(G-E)x'i^)k)v 

-{-G-E)I^J' 

-({GXl^-  GlJ~{-  G-E)XJJk- GI^')u' 

-{{GX^I^-  GXlJ-i-  G-E)x'iJk-  GXJ^  ')u 
-(-Glyy-  Gl^)o" 

-i-GlJ-GlJjo' 

-[El^^-GX^ly,-GX^lJo  =  -(V/«)o>  (4.4) 

The  effects  of  torsion,  curvature,  and  the  area  moment  of  inertia  of  the  cross 
section,  as  well  as  the  variation  of  these  quantities  along  the  rod.  are  next  examined.  When 
these  equations  are  compared  with  the  equilibrium  equations  derived  by  Kingsbury  8.  9 
and  others  for  various  special  cases,  certain  differences  are  observed.  The  effects  of  these 
differences  for  the  case  of  a  circular  ring  au'e  examined  in  Chapter  6. 


18 


4.2  Straight  Rod  with  Variable  Section  Properties 


In  the  case  of  a  straight  rod.  k  and  A  and  their  derivatives  are  equal  to  zero:  thus. 
Equations  (4.1).  (4.2).  (4.3)  and  (4.4)  are  simplified  to  : 


rr  ""  nr<T  >  I"  rr  >'  "  r>T  H"  nor  <  "  " 

El^v  ~2EI^  V  -EI^  V  ~EI^u  u 

=  -  Apu" 


rr  ""  nT?T  I  >>  "  pr  ">>  opt  >  pt  "  " 

EI^v  V  -EI^  V  ~EI^u  -2E1^  u  ~EI^  m 


=  -  Apv' 


-  AEw' A ' Ew'  =  —  Apw" 


(-  Gl-  G/Jo"-(-  GlJ-  GlJ)o'  =  -  (Iy,-lJpo- 


(4.5) 

(4.6) 

(4.7) 

(4.8) 


From  the  above  equations,  it  is  seen  that  if  the  area  moment  of  inertia  (I^^)  is 
equal  to  zero,  then  u  and  v  are  uncoupled. 

4.3  Rod  Curved  in  a  Plane  with  Constant  Curvature  and  Variable  Section 
Properties 

In  this  case,  the  rod  is  initially  in  the  form  of  a  segment  of  a  circle  so  that  /c  is  a 
constant  and  A  is  zero.  The  Euler-Lagrange  equations  are  thus  simplified  to  : 

EI^w"'-2EiJkw’'--{{EIJ’-AE)k-GI^^)w’ 

f  Z  PT  lll>  OPT  f  III  /PT  II  /->T  2\  II  I  2  I 

-  Glyy  k^w^EI^v  ~2EI^v  ~{E1^  -GI^)v  -  GI^  k  v 
PT  III!  OPT  I  III  /PT  II  /-.T  2\  II  /-,T  I  2  I  AO  2 

-Ely^n  -~2E1^u  ~{EI^  -GIy^^)u  -GI^km~AEku 
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~{-G-E)I^o"-{-G-2E)iJko'-EiJ'ko  =  -  Apu  (4.9) 

EI^w"'-2E1^'kw"~{E1^”^-GI^^)w'-G1^'k\ 

-El^v' ' "~2E1Jv  ' GI^^W GjJk\' 
-E1^^""~2E1^'u"^-{E1J'-GJ^-^)u"-G1^'k\' 

G-  E)I^o'  G-2E)J^' Ko'- EIjJ ' KO  =  -  Apv"  (4.10) 

{-EI^^-AE)w''-{-EIJk^-a'e)w'-GI^*w 

j-tr  f  ^  f  n  ^  ^  ^  T  3^ 

-  El^  V  ~  El^  KV  ~  GI^  u 

- El^ u"-EiJku  ' -{Gly^^~AEK)u' -a' Ek u 

-{G~E)Ij^"^o' -El^' K^o  =  -  Apw"  (4.11) 

(-G-  E)I^^w'~GI^'k^w 

-[-G-E\I^v"-GIJkv' 

~{-G-E}I^u"-GJ^'ku' 

-(-  GV  GlJo"-{~  GlJ-  GlJ)o'-EI^^o  =  -  [I^-IJpo-  (4.12) 

The  above  equations  are  fully  coupled  in  the  variables.  However,  if  is  zero  as 
in  the  case  of.  for  example,  a  circular  cross  section,  then  two  pairs  of  equations  result. 
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4,4  Rod  Curved  in  a  Plane  with  Variable  Curvature  but  Constant  Section 
Properties 

The  rod  curved  in  a  plane  with  variable  curvature  but  constant  section  properties 
implies  that  the  A  and  its  derivatives  are  zero  and  that  the  derivatives  of  the  section 
properties  are  zero  also.  Substituting  these  vedues  into  the  general  form  Elquations,  (4.1) 
(4.2)  (4.3)  and  (4.4).  yields  : 

£/y^u:'"-3E/^'w"~(3£/^''-G/^^-AEK)w' 

-{Eiy'-3GI^^K')w 

trr  >>"  2  It  '  ' 

-El^v  -Gl^  V  ~2GI^k  V 

-El^u""-  G1^\"-2GI^k  'u'-AEk^u 

G-  E]I^o" G-2E)I^'o'-  EI^"o  =  -  Apu  (4.13) 

EI^w'  "-3E1^  'w'  '-{3EJ^  ’ GJ^^)w'~{EI^ ' ' 3GI^\' )w 
-EI^v  GI^\' 2GI^k  'v' 

-EI^u''"-G1^\"-2GI^k'm’ 

-[-G-E)I^o"-{-G-2E)I^'o'-E1^"0  =  -Apv  (4.14) 

(-  EI^^-  AE)w' 2EI^k  V^(-  EI^k  ’  '-GI^*)w 
-EI^v"'~GI^^v' 

-  Ely^u"  ~{GIy^^-~AEK)\t  +AEk\ 
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-{G—E)I^‘o'—EJ^k'o  =  -  Apw"  (4.15) 

(-  G-  E)I^^w' -{-2G-  E)I^k'w 
-(-G-  E)l^v"-GI^'v' 

-(-  G-E)I^u"-Gl^'u 

-{-GI^-GIJo"~EI^^o  =  -(V/Jpo-  (4.16) 

If  is  set  equal  to  zero  in  the  above  equations,  the  equations  governing  in-plane 
motions  Eq.  (4.13)  and  Eq.  (4.15)  uncouple  from  those  of  out-of-plane  motions  Eq.  (4.14) 
and  Eq.  (4.16)  . 

4.5  Rod  with  Constant  Curvature  and  Torsion  and  with  Cross  Section 
Varying  along  the  Centroidal  Line  and  Symmetric  to  the  Local  Origin 

In  this  special  case,  the  rod  axis  is  a  helix  and  I^^  is  zero.  The  governing  equations 

become  : 

-{2EI^'-E\I^)kw'' 

-{{ElJ E\I^'~  AE)k- GI^^)w' 

-{GXI^-GlJ)K^rv 

cr  >>" 

-EI^v 

-(-EXI^-2E1J-EXI^)v'” 
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-(-  GI^^-2EXlJ-El^''-EX^I^-EXlJ)v" 
-{(GXI^-GiJ-GXI^)k^-EXiJ'-EX^I^')v' 

-{GXlJ-GX^I^yh 

-El  u'"' 
yy 

-2EJ  'u" 

•'yy 

-i-GI^'^-ElJ'-EXlJ-EXhju'’ 

-{-GiJk^-EXI^''-EX^iJ)u' 

GXI^'-GX^I^-AE)k'^u 

-{-G-E)I^o" 

-(-  GXI^-{-  G-2E)I^'-{G-E)XI„)ko' 

-(-  GXl^'- EI^" -GX^I^-EXiJ)ko  =  -Apu  (4.17) 

rr 

EI^w 

-{EXI^-2EIJ)kw" 

-{{EXI^'-EjJ')k-GI^^)w' 

-{-gxi^-giJ)k\ 

-2ElJv'" 

-(-  GI^^-EX'^I^-EXlJ-ElJ’)v' 
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-{GX^1^-GXI^')k^v 

rr  >"> 

-{EXl^~2ElJ-EXI„)u'" 

-{-GI^^-EXlJ-El^"-BX^I^-2EXlJ)u" 

-((-  GXI^-  GI^'-  GXIJk^-EX^iJ~EXiJ')u' 

-(-GX^1^-GXiJ)k\ 

-{-G-E)I^o" 

-ii-2G-E]Xl^~{-  G-2E)iJ)ko' 

-{-GX^l^-{-G-E)XlJ-ElJ')<o  =  -Apu  (4.18) 

( -  EI^ AEW”-(-E1Jk^- a’ E)w'~ GI^ * w 

-EI^v'" -{EX  I^-  EI^’)k  v"-{  GI^^-EXiJk)v'- GXI^\ 

-EI^u"'-{-EiJ-EXQku" 

-(G1^^-(AE~EXI^')k),i'-(GXI^^-a'Ek)u 

-{G-E)I^^o'-(GXI^-EiJ)k^o  =  -Apw"  (4.19) 

(-G-E)I^^w'-{GXI^-GiJ)k\ 

G-E)l^v"~{{2G-E)XI^-  GiJ)kv’^{GXI^'-  GX^I„)kv 
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(-  GI^-  GIJo" 

(-  GlJ-  GjJ)o'-(EI^^-GX^I^-GX^IJ<i>  =  -  (f^-IJpo-  (4.20) 


CHAPTER  6 


RAYLEIGH-RITZ  SOLUTION  FOR  THE  HELICAL  ROD 

In  this  chapter,  the  general  form  of  the  governing  equations  is  specialized  for  the 
helical  rod  which  has  constant  torsion  and  curvature.  The  strain  energy  function  for  the 
helical  rod  is  then  employed  with  the  Rayleigh-Ritz  method  to  formulate  approximate 
solutions  for  mode  shapes  and  natural  frequencies  both  for  a  helical  rod  with  both  ends 
fixed  (clamped-clamped)  and  for  another  with  one  end  fixed  and  the  other  free 
(clamped-free).  For  the  first  of  these  cases  two  different  sets  of  trial  functions  are  used. 
One  of  these  is  the  set  employed  by  Volterra  and  Morell  17  and  by  Den  Hartog  5  for 
partial  rings  with  clamped  ends,  while  the  second  set  employs  the  normal  mode  functions 
for  clamped-clamped  straight  beams.  Use  of  both  types  of  functions  allows  comparison  of 
frequencies  and  mode  results  obtained  from  the  two  different  sets  of  trial  functions.  This 
comparison  enables  us  to  estimate  the  accuracy  of  the  Rayleigh-Ritz  solution  for  the 
clamped-free  case  which  employs  only  the  normal  mode  functions  of  straight  rods 
(admissible  functions). 
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5.1  Strain  Energy  Function  for  the  Helical  Rod 

The  following  expression  for  the  strain  energy  function  for  a  rod  with  constant 
curvature  and  torsion  is  obtained  from  the  general  strain  energy  expression  derived  in 
Chapter  2. 

Helical  Rod  ~ 

{{E  v"‘^-E u'^-G v'^-G 0^~Ek^  (l>^)  2 
—  G  X  u  V  —  G  K  o'  V  —  E X  u  V  ^ ~  E K  o  V  ^ — E K  o  X  u  —  G K  X  o'  V.) 

—  (Ek  u''  v" -  E X  v  v" —Eu"  v" —Ek  X  w'  u'—  EX^  u'  v' ~EX  u'  u" 

—  E O  w' ~E K  o  X  v' -  E K  o  u"  —  G v'  W~-G  X  U  W~  G  o'  w 

~  G  X  vv'  ~  G  u'  V  —G  K  X  o  V  —  G  X^  U  V—  G  K  X  o'  V 

-Gk^  A  u  u'-Gk  o'  u'-Gk  X^ou) 

-{(Eu"^-EX^  v'^-Ek'^  u/'^-Gk-*  1i;2-G^-* 

-Gk^  A^  v^-Go'^-GA^o^)  2 

~  E  K  w'  X  v' -  G  A  vw—Gk'^  u  w—Gk^  X  o  w—  EX  v'  u" —Ek  w'  u" 

-  G  K^  X  u'  V-  G K  X^  o  V—  G K  X  o  v')  1^ 

-{{Ew'^-Ek^  u^)  2-Ekuw')A  (5.1) 

Although  the  section  properties  in  the  above  strain  energy  function  remain 


functions  of  z.  in  the  subsequent  sections  of  this  chapter  these  section  properties  will  be 
kept  as  constant. 
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5.2  Free  Vibrations  of  a  Clamped-Clamped  Helical  Rod  with  Constant 
Section  Properties 


With  reference  to  Section  3.1,  the  natural  boundary  conditions  for  the 
clamped-clamped  case  are  u  ^  =  0  .  u  =  0  .  v  ^  =  0  ,  v  =  0  ,  w  =  0  and  <t>  =  0  for  z  =  0  ,  z  =  1. 

Following  the  solution  of  \’olterra  and  Morell  and  of  Den  Hartog  for  partied  rings, 
suitable  expressions  for  the  displacements  satisfying  the  above  boundary  conditions  are  : 


u  =  / sin  (flf)  {  fl]  1  -  cos  ( 

-  a2  1  -  cos 

-  Us  1  -  cos  (^^)  } 


(5.2) 


i'  =  Isin  (Df)  {  6]  1  -  cos  ( 


2rr,. 


—  65  1  --  cos  ( 


I 


_  63  1  ~  cos  i^)  } 


w  —  lsm{Ut)  Cj  sin  — C2sin  -  C3  sin  (-  ^ 


,2  -  z, 


o  =sin  (Qt )  dj  sin  (^)  -  ^2  sin  (^-^)  —  sin 


(5.3) 

(5.4) 

(5.5) 


where  a,,  b,.  c,  and  d,  are  coefficients  to  be  determined. 


Upon  substitution  of  these  displacements  into  the  strain  energy  function  (5.1)  and 
the  kinetic  energy  function  (3.4).  an  expression  for  the  total  energy  for  a  rod  undergoing 
free  vibration  is  obtained  : 


(5.6) 
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This  expression  is  not  presented  explicitly  in  this  report  because  of  its  excessive  length. 


The  total  energy  is  then  minimized  with  respect  to  the  unknown  coefficients  by 

requiring 

E.F.  ,  =  0  (5.7) 

where 

i  =  .  Cy  .  C2  .  .  d2  .  d^ 

This  yields  a  s.vstem  of  algebraic  equations  : 

C  =  0  (5.8) 

where 

O  —  <2|  .  fl2  .  <2^  ‘  *  ^2  '  ^3  ■  ■  ^2  ■  ^3  ■  ‘  ^2  ^3 

E 

where  H  is  the  natural  frequency  of  the  free  vibration  of  the  rod. 


M(j.  Mj,  and  the  integrated  functions,  col  to  co58.  which  appear  in  the  matrix,  are 
shown  in  Appendix  B  in  FORTRAN  codes.  Since  the  coefficients  of  the  column  matrix  C 
cannot  be  trivial,  characteristic  equations  for  the  clamped-clamped  case  are  obtained  from 
the  condition  : 

,  M„- =0  (5.9) 

Next,  the  same  process  is  repeated  using  normal  modes  of  clamped-clamped 


29 


beams  and  bars  as  comparison  functions. 


As  obtained  from  Young  and  Felger  20  ,  the  normal  mode  functions  are  : 


u  =  /  sin  (fit)  {  fli  {  cosh(^^)  -  cos  (■^^)  ~  kkl  sinh('^^)  -  sin  } 

-  02  {  cosh(-^^)  -  cos  (•^^)  -  kk"!  sinh(-^^)  -  sin  }}  (5.10) 


V  =  /sin  (Jit)  {  ij  {  cosh(^^)  -  cos  -  kk\  [ sinh(.^^)  -  sin  (■^^)  } 

—  62  {  cosh(-^^)  -  cos  -  1:1:2  sinh(.^^)  -  sin  }}  (5-11) 


(5.12) 


o  =  sin  (fl  t)  </j  sin  (—)  —  sin  (^-~) 


(5.13) 


The  coefficients  of  kl.  kkl.  k2  and  kk2  3  .  which  appear  in  Equation  (5.10) 


through  (5.13)  are  given  in  Appendix  C.  Upon  substituting  Equation  (5.10)  through 


Equation  (5.13)  into  the  kinetic  and  strain  energy  functions,  and  repeating  the 


minimization  process,  another  set  of  characteristic  equations  is  obtained. 


M,  -  $  =  0 


(5.14) 


where  the  elements  of  the  matrices  and  the  coefficients  are  given  in  Appendix  C. 


Frequency  results  for  an  ‘‘edmost  straight"  helical  rod  for  the  two  sets  of  trial 


functions  are  compared  with  each  other  and  with  straight  beam  natural  frequencies  in 


Table  5-1. 
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Helical  Rod 
Eqs.  5-2. 5-3 

5-4.  and  5-5 

Helical  Rod 

Normal 

Mods  Functions 

straight  Beam 

u 

0.226*7  X  10^ 

0.2263  X 

104 

0.2263 

X  10^ 

Ist 

V 

0.2267  X  10^ 

0.2263  X 

104 

0.2263 

X  10^ 

Mode 

w 

0.6355  X  10® 

0.6355  X 

lo® 

0.6355 

X  10® 

♦ 

0.4019  X  10® 

0.4019  X 

10® 

0.4019 

X  10® 

u 

0.1234  X  105 

0.6238  X 

10< 

0.6238 

X  10^ 

2nd 

V 

0.1234  X  10^ 

0.6238  X 

10* 

0.6238 

X  10^ 

Mode 

w 

0.1271  X  lO”^ 

0.1271  X 

lo"^ 

0.1271 

X  lo"^ 

0.8038  X  10® 

0.8038  X 

10® 

0.8038 

X  10® 

Unit 

.  of  Freq.“  rad/aec.  Length  of 

The  Rod 

■  1  inc 

G/E  -  0.4,  LK  -  1x10“®,  LL  - 
Density  =  0.283  Ib/inch^ 

1x10* 

■®,  LD  - 

lxl03. 

Table  6-1:  Frequencies  of  Clamped-Clamped  Straight-like  Rod 

For  this  structure,  the  first  modes  of  the  natural  frequencies  are  almost  identical 
for  all  three  set  of  results.  The  fi'equencies  corresponding  to  the  second  set  show  that  the 
beam  normal  mode  function  approximation  for  the  helical  rod  with  small  curvature  gives 
frequencies  more  nearly  equal  to  those  of  the  straight  beam  and  are  therefore  more 
accurate.  Because  the  hyperbolic  ^e  and  come  functions  are  mode  functions  for  the 
strught  beam,  their  use  will  give  poor  results  for  rods  with  large  curvature.  Thus,  the  first 
set  of  admissible  functions  would  be  expected  to  give  better  results  for  natural  fiequencies. 
That  this  is  the  case  will  be  shown  in  Chapter  6. 
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5.3  Free  Vibrations  of  a  Clamped-Free  Helical  Rod  with  Constant  Section 
Properties 


In  the  clamped-free  case,  the  boundary  conditions  are  =0.  u  =  0,  =  0,  v  =  0. 

w  =  0.  o  =  0  at  2  =  0  and 


dF 

a  u 


0 


d  dF  dF 


d  u 


d  zd  u 
dF 


77 


=  0 


T7 


=  0 


dv 

^  dF  _  dF 
d  zd  V  ’  d  V 

dF 


=  0 


d  If 

dF 

Fo 


=  0 


=  0 


zX  z  -  1. 


The  relationships  between  the  above  expressions  and  the  force  and  moment 
resultants  acting  on  a  cross  section  are  not  obvious.  In  Kingsbury's  report  8.9  .  the 
moments  and  forces  are  shown  to  be  related  to  the  displacement  components  in  the 
following  manner  : 


M,=  -  El^aJ  -KO) 

(5.15) 

My=El^ay  -  El^iaJ  -K<t>) 

(5.16) 

A/,=  G/aXOj- G/^xOj, 

(5.17) 
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=  Ku)  (5.18) 

V^=-M^'-XM,  (5.19) 

V^=mJ  -  XMy-KM,  (5.20) 

where 

I  , 

a^  =  u  -  X  V  —  K  w 
I  > 

Oy  =  -  t)  -  A  ti 

By  substituting  and  into  Equation  (5.15)  through  Equation  (5.20),  the 


following  equations  are  obtained  : 

M^=-  EI^{u"  -  Xv'  ~Kw' )~  EI„{-v"  -  Xu'  -Ko)  (5.21) 

My  =  E  l^{u"  -  X  V  -  K  w' )  -  Eljy{-v"  -  X  u  -  K  o)  (5.22) 

M,=  G{I^-l^]o'  -G1^k{v'  -Xu)-  G1^k(u'  -Xv-^w)  (5.23) 

V^=EA{w'  -  Ku)  (5.24) 

V^=-  EIyy(u  -Xv  -ku)  )-EI^(-v  -Xu  -ko  ) 

—  X  E  I^{u"  -  X  v'  —  K  ui' )  -  X  E  -  v"  -  Xu'  —  Ko)  (5.25) 

Vy  =  -  EI^{u  -Xv  —KW  )  —  Eljg{-v  -Xu  -Ktp  ) 

-  X  Blyy  {u"  -  X  v'  —  K  w'  )  -r  X  E  I^{- v"  -  X  u'  K(l>) 


-GK{1„-Iyy)<t>'  -GI„k^{v’^Xu)^GI^k^u'-  Xv  +  Kw)  (5.26) 


Upon  substitution  of  the  strain  energy  function,  Equation  (5.1).  boundary 


conditions  at  the  free  end  are  found  to  be  expressed  in  terms  of  moment  eind  force 
resultants  in  the  following  manner: 


|^  =  A/,=  0 

(5.27) 

o  r 

(5.28) 

a  w 

(5.29) 

o 

II 

II 

CO  |(^ 

(5.30) 

d  zd  u  d  u 

(5.31) 

d  dF  _ 

d  zd  v"  d  V  ^ 

(5.32) 

It  is  concluded  that  the  forces  and  moments  are  zero  oru  =  u^  =  u^^  =  \i'"  =  v 
=  \'  =  v"  =  \-" '  =  w  =  w'  =  w"  =  o  =  =  0  at  the  free  end  [z  =  /).  In  this  situation, 

it  is  very  hard  to  find  the  suitable  functions  which  fully  satisfy  the  boundary  conditions. 
Hence,  the  normal  mode  functions  of  a  clamped-free  straight  beam  which  are  admissible 
trial  functions  are  employed  to  get  approximate  solutions. 


The  normal  mode  functions  for  the  straight  clamped-free  beam  are  : 


u  =  / sin  (fit)  {  Uj  {  cosh(^^)  -  cos  -  kkl  [  sinh(-^^)  -  sin  (■^^)  ]} 

—  02  {  cosh(^^)  -  cos  “  kk2  '  sinh(^^)  -  sin  (■^^)  ]}}  (5.33) 


When  the  torsion  \  and  curvature  k  become  small,  the  helical  rod  will  approach  a  straight  beam  and 
the  boundary  conditions  will  be  much  the  same  as  those  of  the  straight  beam  (u**  =  a' "  =  v"  =  = 


w  '  =  o'  =  0 ) . 
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v=  Ism  (n<)  {  {  cosh(^^)  -  cos  -  kkl  sinh(^^)  -  sin  } 

—  62  i  cosh(-^^)  -  cos  -  kk2  sinh(^^)  -  sin  (■^^)  }}  (5.34) 


w 


=  I  sin  (n  t)  Cl  sin  (^)  -  Cj  sin  (^) ' 


(5.35) 


o  =  sin  (n  t)  '  di  sin  (|^)  -  (ij  sin  (^)  : 


(5.36) 


By  substituting  the  above  expressions  into  the  potential  and  kinetic  energy- 


functions.  and  repeating  the  steps  as  in  the  last  section,  the  following  equation  is  obtained: 


M,  -  $  Mf  =  0 


(5.37) 


where,  the  elements  of  the  .M^.  matrices  and  the  coefficients  are  shown  in  Appendix  D. 


The  frequencies  for  clamped-free  helical  rods  with  circular  cross  sections,  small 


curvature,  and  small  torsion  are  tabulated  in  Table  5-2.  These  results  are  compared  and 
found  consistent  with  those  of  straight  rods. 


Hellral  Rod 

Straight  ] 

Normal 

Mode  Functions 

u 

0.3556 

X 

10® 

0.3556  X 

10® 

1st 

V 

0.3556 

X 

10® 

0.3556-X 

10® 

Mode 

w 

0.3177 

X 

10® 

0.3177  X 

10® 

♦ 

0.2010 

X 

10® 

0.2010  X 

10® 

u 

0.2229 

X 

10® 

0.2229  X 

10® 

2nd 

V 

0.2229 

X 

10® 

0.2229  X 

10® 

Mode 

V 

0.9532 

X 

10® 

0.9532  X 

10® 

4* 

0.6029 

X 

10® 

0.6029  X 

10® 

Unit  of  Freq.“  rad/sec.  Length  of  The  Rod  ■  1  inch 
G/E  -  0.4,  LK  -  1x10"®,  LL  -  1x10"®,  LD  -  1x10®, 
Density  *  0.283  Ib/inch® 


Table  6-2:  Frequencies  of  Clamped-Free  Straight-like  Rod 


CHAPTER  6 


COMPARISON  OF  VARIOUS  EQUATIONS  AND  THE 
STRAIN  ENERGY  FUNCTIONS  FOR  PARTIAL 
RING  WITH  CIRCULAR  CROSS  SECTION 

Through  comparison  of  the  equations  of  motion  of  curved  rods  derived  in  this 
report  vi-ith  those  derived  earlier  by  Kingsbury  8.  9  .  slight  differences  are  noted  between 
the  two.  These  differences  do  not  represent  errors  in  one  or  the  other  formulation  but 
result  instead  from  the  way  in  which  the  equilibrium  and  kinematic  terms  are  combined 
and  the  order  in  which  the  various  simplifying  assumptions  are  invoked.  This  chapter 
assesses  the  effect  of  these  differences  on  the  vibration  of  a  circular  ring  which  is  a  special 
case  of  general  space  curvature  for  which  the  torsion  is  zero  and  the  curvature  is  constant. 

The  present  equations  are  first  compared  with  those  of  Kingsbury  and  with  those 
presented  by  A.  E.  H.  Love  10  for  the  circular  ring.  It  is  found  that  the  most  essential 
differences  appear  in  the  equations  governing  out-of-plane  vibrations  with  those  of 
Kingsbury  and  Love  being  most  nearly  in  agreement. 

The  quantitative  effect  of  these  differences  is  assessed  by  employing  the  suitably 
simplified  version  of  the  strain  energy  functions  derived  in  chapter  2  with  a  corresponding 
one  presented  by  Love  to  obtun  Rayleigh-Ritz  solutions  for  the  natural  frequencies  of 
out-of-plane  vibrations  of  a  clamped-clamped  ring. 
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6.1  Comparison  of  Force  Equilibrium  Equation 


6.1.1  Derivation  of  Governing  Equations  using  Strain  Energy  Function 


Upon  applying  the  equations  given  in  Section  4.5,  eliminating  all  the  terms 
containing  torsion  A.  or  their  derivatives,  and  dropping  the  terms  containing  the 
derivatives  of  the  area  moment  of  inertias  and  area  A,  the  following  equations  are 


obtained  : 

-  iAEK-Gly^^)w' -  Gly^^u" -AEk^u  =  -  Apu  (6.1) 
-{EI^^-AE)u-"~GIy^*w-EI^J”-{GI^^-AEK)u’  =  -  Apw  (6.2) 

EI^v""-Gl^'^v'-[G-E)I^o"  =  -  Apv  (6.3) 

-{G-E)l^v"-iI^-IJGo”~EI^^o  =  -  (V/Jpo-  (6.4) 


It  is  noted  that  for  the  ring  problem,  the  equations  governing  in-plane  motions 
(Eq.  (6.1)  and  Eq.  (6.2))  uncouple  from  those  of  out-of-plane  motions. 

6.1.2  Governing  Equations  from  Kingsbury’s  Formulation 


EI^w"'-AEKw'-EIyyu""-AEK^u  =  -  Apu  (6.5) 

-{EI^'^-AE)w"-EI^u"-AEku  =  -  Apw"  (6.6) 

ElJ'"-Giyv''-[G{I„-^I^)-EI„)K<i>''  =  -  Apv  (6.7) 

-{G~E)I^v"- [I^-IJGo" ^EI^'^6  =  -  {Iyy--Ij4>"p  (6.8) 
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Through  comparisons  of  Equations  (6.1)  and  (6.2)  with  (6.5)  and  (6.6),  differences 
between  these  two  pairs  are  seen  to  involve  terms  of  the  order  of  the  squaire  or  higher  power 
of  the  curvature.  These  terms  should  not  have  a  significant  effect  on  the  results  for  rods 
with  small  curvature.  Equation  (6.7)  contains  an  extra  term  in  the  coefficient  of  which 
causes  its  magnitude  to  be  approximately  twice  that  of  the  corresponding  coefficient  of 
Equation  (6.3).  The  practical  effect  of  this  difference  can  only  be  assessed  numerically. 
Equations  (6.3)  and  (6.4)  lead  to  a  symmetric  stiffness  matrix  in  the  characteristic  equation 
for  the  eigenvalues  while  Equations  (6.7)  and  (6.8)  do  not. 

6.1.3  Governing  Equations  based  on  Love’s  Assumptions 

By  using  Love's  equations  of  motion,  stress  resultants,  and  couples  of  the 
3 

ring  10  .  the  out-of-plane  vibration  equations  can  be  formulated  as  shown  below  : 

EI^v""-GiI^-ly^)K\"-{G{l^~IJ~EIjKo"  =  -  Apv-  (6.9) 

-iG~E)I^v"-(I^-IJGo"-EI^^o  =  -  {ly,~IJpo-  (6.10) 

Equations  (6.8)  and  (6.10)  seem  to  be  identical,  while  Equations  (6.7)  and  (6.9) 
differ  only  by  the  difference  between  and  It  may  be  concluded  that  Love's  equations 
for  out-of-plane  vibration  more  clearly  resemble  those  of  Kingsbury. 

In  order  to  compare  Equations  (6.1)  and  (6.2)  with  the  corresponding  in-plane 
vibration  equation  derived  by  Love,  it  is  necessary  to  simplify  the  former  equations  by 
imposing  the  condition  of  inextensibility  of  the  ring  axis.  The  vanishing  of  axial  strain 

3 

Equations  18.  19  and  20.  p.  451. 


implies  : 
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d  w 
~dz 


=  ^■  li 


Upon  substituting  the  above  relation  into  Equations  (6.1)  and  (6.2),  deleting  the 
shear  terms,  then  changing  the  independent  variable  from  z  to  B,  which  is  the  angle  of  the 
arc  of  the  ring,  the  following  equations  are  obtained  : 


yy 


^6 

aO  W 

K  - X  - 

- 

4  d*  w 

h.  . 

dPT 

yy 

dB* 

4^^  Ml 

t\.  . 

48^  w 

K  - 5 

dB' 

yy 

dB^ 

i=- pA 


d*  w 

d¥J? 


=  pA 


w 


(6.11) 

(6.12) 


Adding  Equations  (6.11)  and  (6.12)  yields  the  same  equation  as  given  in  Love's 
book  (page  452).  If  we  integrate  Equation  (6.11)  twice  with  respect  to  the  variable  6.  then 
add  the  result  to  Equation  (6.12).  we  also  find  the  constraints  for  stress  resultants  and  the 
flexural  couple  vanish  at  the  ends  when  the  ring  is  incomplete^  10  .  These  results  show 
that  the  equation  found  by  Love  can  be  obteiined  by  simplifying  and  combining  in-plane 
vibration  Equations  (6.1)  and  (6.2)  through  the  relation  of  inextension  of  the  ring  axis. 


4 

p  452.  “When  the  ring  is  incomplete  the  frequency  equation  is  to  be  obtained  by  forming  the  conditions 

^  ^  I  ■  ^  ^ 


that  N.  T.  G  vanish  at  the  ends” 


6.2  Comparison  of  Strain  Energy  Functions 


6.2.1  Out-of-Plane  Vibration  Strain  Energy  Function  of  Circular  Ring 

To  obtain  the  strain  energy  for  the  circular  ring  from  Equation  (2.17),  A,  and 
the  derivatives  of  k  are  set  equal  to  zero.  This  results  in  the  following  expression  : 

^■^■^■CiTcularRtng  Chd-  of-  PlaneVibrction  ~ 

|4(r")2  -  EI^k  o  f"  (v')^ 

-  GI^ko' v' (6.13) 

6.2.2  Out-of-Plane  Vibration  Strain  Energy  Function  of  Circular  Ring 
Derived  from  Love’s  Formulation 


.\ccording  to  Love  10  .  the  total  strain  energy  of  the  circular  ring  is  : 


dz) 


(6.14) 


EI^  G(/„-/„) 

where  G  and  H  are  the  torsional  and  bending  moment®  10  given  by  Equations  (6.15)  and 
(6.16). 


G  —  E o  -  V 


H=G(;„-/^)  {K  v'~o') 


(6.15) 

(6.16) 


Substitution  of  Equations  (6.15)  and  (6.16)  into  equation  (6.14)  yields  the  strain 


equation  20.  p.4ol . 
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energy  function  (6.17). 

S-E.F.  CircuiarRmg  Chii—  of—  Plant  Vibratton  ~ 

V  -  o  v"  -  i  G  (/„ -  1^)  (^')2 

(6.17) 

Comparing  the  two  strain  energy  functions.  Equations  (6.13)  and  (6.17).  w’e 
conclude  that  the  coefficients  of  the  (v*)^  and  o'v'  terms  differ  by  a  factor  of  two  if  = 
lyy.  The  quantitative  effect  of  these  differences  on  the  lowest  frequency,  f.  for  the 
out-of-plane  vibration  of  the  clamped-clamped  ring  is  tabulated  in  Table  6-1. 

In  Table  6-1.  the  first  and  second  sets  of  data  are  the  results  of  the  substitution  of 
admissible  functions,  given  by  Equation  (5.2)  through  Equation  (5.5).  into  the  strain  energy 
expressions  (6.13)  and  (6.17)  respectively.  For  the  purpose  of  comparison,  the  results  of  the 
substitution  of  the  normal  mode  functions  (admissible  functions)  into  Equation  (6.17)  are 
shown  in  the  third  set  of  data. 

The  three  sets  of  data  will  be  compared.  In  the  first  set  of  data,  results  for  6  > 
180°  are  the  same  as  the  data  presented  by  Den  Hartog  [o/  .  Comparing  the  first  and 
second  sets  of  data,  it  is  found  that  the  differences  of  the  two  sheeu"  strain  energy  terms  will 
affect  /  to  the  maximum  15  %.  Here,  we  may  conclude  that  the  two  formulations  can  yield 
different  results  for  large  included  angles,  but  there  is  no  way  of  assessing  their  relative 

V(2,''(DCfl))V{*) 


'Oen  Hartog's  energy  function  is  the  same  as  equation  (6.17) 
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Equation 

(6-17) 

Equation 

(6-13) 

Equation 

(6-13) 

Admissible  functions 

Normal 

Mode 

Eqs.  5- 

2,  5-3,  5- 

•4  and  5-5 

Functions 

DD-0.02 

DD-0.1 

DD-0.02 

DD-0.1 

DD-0.02 

DD-0.1 

e 

f 

f 

f 

f 

f 

f 

10 

734.2 

227.8 

734.8 

227.9 

733.7 

227.9 

20 

183.0 

113.9 

183.3 

114.2 

183.2 

114.2 

30 

80.87 

73.12 

81.19 

74.36 

81.19 

74.36 

40 

45.14 

44.06 

45.45 

44.75 

45.51 

44.81 

50 

28.62 

28.18 

28.92 

— 28.64 

29.00 

28.73 

60 

19.66 

19.42 

19.95 

19.80 

20.04 

19.90 

70 

14.27 

14.12 

14.54 

14.45 

14.65 

14.56 

eo 

10.78 

10.69 

11.05 

10.99 

11.16 

11.10 

90 

8.410 

8.337 

8.654 

8.613 

8.731 

100 

6.708 

6.662 

6.949 

6.920 

7.041 

110 

5.464 

5.431 

5.693 

5.672 

5.795 

120 

4.524 

4.500 

4.742 

4.726 

4.852 

130 

3.799 

3.781 

4.006 

3.994 

4.121 

140 

3.230 

3.216 

3.425 

3.416 

3.544 

150 

2.774 

2.764 

2.959 

2.954 

3.081 

160 

2.406 

2.398 

2.581 

2.575 

2.705 

no 

2.104 

2.097 

2.269 

2.265 

2.252 

180 

1.853 

1.848 

2.009 

2.006 

2.137 

190 

1.644 

1.64  0 

1 .791 

1.789 

200 

1.4  67 

1.463 

1.606 

1.604 

210 

1.316 

1.314 

1.448 

1.447 

220 

1.187 

1.185 

1.312 

1.311 

230 

1.076 

1.074 

1.195 

1.194 

240 

.97  95 

.9780 

1.092 

1.091 

250 

.8954 

.8941 

1.002 

1.001 

260 

.8216 

.8206 

.9231 

.9224 

270 

.7567 

.7559 

.8531 

.8526 

280 

.6994 

.6987 

.7910 

.7905 

X 

e  / 

290 

.6486 

.6479 

.7356 

.7352 

X 

.  / 

300 

.6034 

.6026 

.6861 

.6857 

\/ 

310 

.5631 

.5626 

.6416 

.6413 

320 

.5270 

.5266 

.6016 

,6013 

330 

.4948 

.4944 

.5654 

.5652 

34  0 

.4659 

.4655 

.5327 

.5325 

350 

.4399 

.4396 

.5301 

.5029 

360 

.4166 

.4163 

.4761 

,4760 

DD  ■ 

>  Diameter 

of  Cross 

Section  / 

Diameter 

of  Ring 

Tabic  6-1;  Frequencies  of  Out-of-Plane  Vibration, 
Clamped  Two  Ends  Partial  Ring  (G/E=0.4) 
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accuracy  from  the  published  literature.  From  the  other  viewpoint,  it  is  found  that  when 
the  angle  of  the  ring  is  small,  the  error  of  the  coefficient  /  between  the  first  and  second  sets 
of  data  is  small  also.  From  the  factor  of  two  mentioned,  it  may  be  concluded  that  when  the 
angle  of  the  ring  is  small,  the  shear  strain  terms  do  not  overly  affect  the  accuracy  of 
frequency. 

From  Table  6-1,  we  also  found  that  most  of  the  frequency  parameter  (/)  values 
obtained  by  using  normal  mode  functions  of  straight  rods  are  slightly  higher  than  those 
which  were  obtained  from  the  other  set  of  admissible  functions.  Missing  data  are  negative 
eigenvalues  which  probably  result  from  use  of  straight  beeim  mode  shapes  and  from  the 
truncation  errors  of  the  IMSL  program. 
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EIGZF  subroutine.  IMSL  INC,  Huston.  Texas. 


CHAPTER  7 


VIBRATION  OF  THE  CLAMPED-FREE  HELICAL  RODS 

In  this  chapter,  the  Rayleigh-Ritz  formulation  of  Chapter  5  for  clamped-free 
helical  rods  is  used  to  examine  the  effects  of  curvature,  torsion,  cross  section  shape  and 
orientation  of  cross  section  principal  axes  on  the  fundamental  natural  frequencies  eind  mode 
shapes  of  a  helical  rod  of  elliptical  cross  section.  In  all  cases  the  rod  material  is  assumed  to 
have  a  Poisson's  ratio  0.3  (GE=G  E=0.4). 

7.1  Effect  of  Curvature  and  Torsion  on  Mode  Shapes  and  Natural 
Frequencies  for  Rod  with  Uniform  Circular  Cross  Section 

In  this  section,  the  effects  of  slenderness,  curvature,  and  torsion  on  the  lowest 
frequency  and  mode  shape  are  examined  on  rods  of  constant  circular  cross  section.  With  a 
fixed  value  of  LD  {I  radius  of  cross  section).  LK  [I  k)  and  LL  (/  >^A)  are  varied.  LD  is 
then  again  fixed  at  a  new  value  and  the  calculations  repeated.  The  results  of  these 
calculations  for  the  fundamental  natural  frequency  of  the  fixed-free  helix  are  presented  in 
Figures  7-1  and  7-2. 

As  expected,  the  fundamental  frequency  increases  with  decreasing  LD.  It  is  seen 
that  the  natural  frequency  is  insensitive  to  changes  in  the  curvature  and  torsion  parameters 
for  small  values  of  these  parameters.  This  effect  is  explored  in  more  detail  in  Figure  7-2 
which  plots  frequency  in  a  linear  rather  than  logarithmic  scale.  It  is  seen  that  the 
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fundamental  natural  frequency  becomes  noticeably  dependent  upon  curvature  and  torsion 
when  each  of  these  parameters  reaches  the  value  of  10“ 

The  effects  of  changes  of  curvature  and  torsion  on  mode  shape  are  explored  in 
Table  7-1.  At  small  v£ilues  of  both  torsion  and  curvature,  transverse  displacement 
dominates  the  mode  shape.  As  the  curvature  is  increased,  for  constant  small  torsion,  the 
displacements  u  and  w  become  coupled  as  do  r  and  o.  On  the  other  hand,  as  the  torsion 
parameter  is  increased  with  small  curvature  (moving  down  the  left  hand  column)  the 
lateral  displacements  u  and  v  become  coupled  and  approach  the  same  magnitude.  If  both 
parameters  are  increased,  then  all  of  the  variables  approach  the  same  magnitude. 
Significant  coupling  among  the  displacement  variables  seems  to  begin  to  occur  at  about  the 
same  values  of  the  curvature  and  torsion  pareimeters  when  frequency  changes  become 
noticeable. 

7.2  Effect  of  Non-Circular  Cross  Section  on  Mode  Shapes  and  Frequencies 

In  order  to  examine  the  effect  of  out-of-roundness  on  curved  rod  motion,  the  cross 
section  of  the  rod  is  changed  from  circular  to  elliptical  while  keeping  the  cross  section  area 
constant.  The  configuration  is  shown  in  Figure  7-3.  The  product  of  inertia  (I^)  is  then 
systematically  varied  by  rotating  the  principal  axes  of  the  cross  configurations  while  the 
torsion  (LL)  and  curvature  (LK)  parameters  are  varied  for  each  value  of  orientation  of  the 
principal  axes  of  the  cross  section. 

Figures  7-5  £md  7-6  show  the  effect  on  fundamental  natural  frequency  of  changing 
the  ratio  of  major  to  minor  axes  while  keeping  .3  =  0  (1^=0).  Figure  7-4  shows  the  results 
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Figure  T>1:  Frequencies  of  Clamped~Free  Helical  Rod. 
{LD=lxl0‘  toLD=lxiO*) 
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Table  7*1:  Fundamental  Frequencies  and  Mode  Shapes  of  Clamped-Free 
Helical  Rod.  (LD=lx  10*,  a/b=l) 
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for  a  circular  cross  section  to  provide  a  basis  for  comparison. 

It  is  seen  that  the  effect  of  out-of-roundness  is  not  synametrical  because  the 
direction  of  the  major  axis  substantially  affects  the  manner  in  which  the  fundamental 
natural  frequency  is  affected  by  changes  of  the  torsion  and  curvature  parameters. 

Displacement  coupling  caused  by  the  non-axisymmetric  cross  section  is  explored 
in  Tables  7-2  through  7-4  and  Figures  7-5.  7-6.  and  7-7.  Figures  7-5,  7-6.  and  7-7  illustrate 
the  effects  of  cross  section  rotation  on  fundamental  naturzJ  frequency,  .\gain.  it  is  noted 
that  rotating  the  cross  section  of  fixed  geometry  changes  the  manner  in  which  the  curvature 
and  torsion  parameters  affect  the  frequency. 

Comparing  Tables  7-2  through  7-4.  it  is  first  noted  that  for  the  almost  straight  rod 
(LL=LK=10~®)  the  dominant  displacement  remains  in  the  direction  of  the  miner  axis  as 
the  cross  section  is  rotated.  For  higher  values  of  curvature  and  torsion,  rotation  of  the  cross 
section  is  seen  to  change  the  relative  values  and  phase  of  the  four  displacement 
components. 

Finally,  it  is  noted  that  when  LD=1000,  LK  and  LL  are  smaller  than  10“^  and 
a  b  is  in  the  range  1.1  to  1.  then  the  fundamental  frequencies  are  independent  of  the 
rotation  angle  3  and  also  lower  than  the  frequencies  of  a  rod  with  circular  cross  section. 


Figure  7*8:  Defonnation  of  Cross  Section 


Figure  7-4:  Effects  of  Non-Circular  Cross  Section  On  the  Fundamental 

Frequency  {a/b=l) 


Figure  7*6:  Ejects  of  Non-Circular  Cron  Section  On  the  Fundamental 
Frequency  (a/b=l.l,^=0’ ) 
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16160 

16160 

16190 

18340 

39170 

U 

O.lOE-fOl 

O.lOE-l-01 

0.10E4-01 

O.lOE-t-01 

O.lOE-t-01 

O.lOE+01 

V 

-0.23E-02 

-0.23E-02 

-0.23E-02 

-0.26B-02 

-0.5eE-t-02 

0.69E+02 

N 

0.54E-08 

0.54E-06 

0.54E-04 

0.54E-02 

0.91E-^00 

0.94E-f02 

0 

0.43E-10 

0.43E-08 

0.43E-06 

0.40E-04 

-0.46E-f02 

-O.13E<t-02 

Table  7-2:  Effects  of  Non-Circular  Cross  Section  On  the  Fundamental 
Frequency  and  Mode  Shape  (a/b=l.l,  /9=0* ) 
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LL  O.lOE-07 


LK 

O.lOE-07 

O.lOE-05 

O.lOE-03 

O.lOE-01 

O.lOE+01 

O.lOE+03 

FREQ 

339.1 

339.1 

339.2 

356.2 

364.0 

10730 

U 

O.lOE+01 

O.lOE+01 

O.lOE+01 

O.lOE-t-01 

O.lOE+01 

O.lOE-t-01 

V 

O.lOE+01 

O.lOE+01 

O.lOE+01 

0.94E<t-02 

0.31E+05  -0.17E+06 

W 

0.39E-08 

0.39E-06 

0.39E-04 

0.54E-02 

0.90E-t-00 

-0.23E+03 

0 

0.38E-08 

0.36E-06 

0.38E-04 

0.39E-«-00 

O.lOE-t'05 

O.74E<»-04 

LL 

O.lOE-05 

LK 

O.lOE-07 

O.lOE-05 

O.lOE-03 

O.lOE-01 

0.10E401 

0.10E-«-03 

FREQ 

339.1 

339.1 

339.2 

356.2 

364.0 

10730 

U 

O.lOE+01 

O.lOE+01 

O.lOE-fOl 

O.lOE-fOl 

O.lOE-t-01 

O.IOE+Ol 

V 

O.lOE-t-01 

O.lOE+01 

0.10E4-01 

0.94E<»-02 

0.31E-t-05 

-0.17E+06 

W 

0.39E-08 

0.39E-06 

0.39E-04 

0.54E-02 

0.90E-«-00 

-0.23E+03 

0 

0.38E-08 

0.38E-06 

0.38E-04 

0.39E-»^00 

0.10E405 

O.74E+04 

LL 

O.lOE-03 

LK 

O.lOE-07 

O.lOE-05 

O.lOE-03 

O.lOE-01 

O.lOE-l-01 

O.lOE-t-03 

FREQ 

339.1 

339.1 

339.2 

356.2 

364.0 

10730 

U 

O.lOE+01 

O.lOE+01 

O.lOE+01 

O.lOE+01 

O.lOE-t-01 

O.lOE+01 

V 

O.lOE+01 

O.lOE+01 

O.lOE+01 

0.95E+02 

0.31E+05 

-0.17E+06 

w 

0.39E-08 

0.39E-06 

0.39E-04 

0.54E-02 

0.90E-)-00 

-0.23E+03 

0 

0.38E-08 

0.38E-06 

0.38E-04 

0.39E-«-00 

O.lOE+05 

O.74E+04 

LL 

O.lOE-01 

LK 

O.lOE-07 

O.lOE-05 

O.lOE-03 

O.lOE-01 

O.lOE+01 

O.lOE+03 

FREQ 

339.1 

339.1 

339.2 

356.4 

364.1 

10730 

U 

0,10E+01 

O.lOE+01 

O.lOE+01 

O.lOE+01 

O.lOE+01 

O.lOE+01 

V 

0.99E+00 

0.99E+00 

O.lOE+01 

O.lOE-t-03 

0.39E<t-05 

-0.22E+06 

W 

0.39E-08 

0.39E-06 

0.39E-04 

0.51E-02 

0.89E-t-00 

-0.32E+03 

0 

0.38E-08 

0.38E-06 

0.38E-04 

0.43E-)-00 

0.13E+05 

O.94E-f04 

LL 

O.lOE+01 

LK 

O.lOE-07 

O.lOE-05 

O.lOE-03 

O.lOE-01 

O.lOE+01 

O.lOE+03 

FREQ 

366.9 

366.9 

367.1 

395.6 

436.2 

10850 

U 

O.lOE+01 

O.lOE+01 

O.lOE-t-Ol 

O.lOE-i-01 

0.10E<f01 

O.lOE+01 

V 

0.44E-t-00 

0.44E-t-00 

0.45E+00 

-0.14E4-02 

-0.17E+04 

0.85E4-04 

W 

0.38E-06 

0.36E-06 

0.38E-04 

0.94E-02 

0.91E<«-00 

0.13E+03 

0 

0.19E-08 

0.19E-06 

0.19E-04 

-0.78E-01 

-0.67E+03 

-O.35E+03 

LL 

O.lOE+03 

LK 

O.lOE-07 

O.lOE-05 

0 . lOE-03 

O.lOE-01 

O.lOE-l-01 

O.lOE-i-03 

FREQ 

>  16160 

16160 

16160 

16170 

17540 

38960 

U 

O.lOE-fOl 

O.lOE+01 

O.lOE+01 

O.lOE-t-01 

O.lOEi-01 

0.10E4-01 

V 

-O.lOE+01 

-O.lOE+01 

-O.lOE+01 

-0.10E-«>01 

-0.25E'l-02 

0.69E<f02 

W 

0.54E-08 

0.54E-06 

0.54E-04 

0.54E-02 

O.BBS+OO 

0.94E4-02 

0 

-0,65E-08 

-0.65E-06 

-0.65E-04 

-C.66E-02 

-0.1BE4-02 

O.12E-I-02 

Tabic  T-8;  Effects  of  Non-Circular  Cross  Section  On  the  Fundamental 
Frequency  and  Mode  Shape  (a/b=l.l,  ^=45* ) 
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LL 

O.lOE-07 

LK 

O.lOE-07 

O.lOE-05 

O.lOE-03 

O.lOE-01 

O.lOE-t-01 

O.lOE+03 

FREQ 

339.1 

339.1 

339.3 

373.0 

379.0 

11260 

U 

O.lOE+01 

O.lOE-fOl 

O.lOE+01 

O.IOE+Ol 

O.lOE-t-Ol 

O.IOE+Ol 

V 

•0.47E-07 

-0.47E-07 

-0.47E-07 

-0.51E+09 

-O.llE+12 

0.59E-»-12 

W 

0.39E-0B 

0.39E-06 

0.39E-04 

0.71E-02 

0.90E-»-00 

0.96E4-01 

0 

0.56E-16 

-0.17E-13 

-0.17E-11 

-O.23E+07 

-0.40E+11 

-0.25E+11 

LL 

O.lOE-05 

LK 

O.lOE-07 

O.lOE-05 

O.lOE-03 

O.lOE-01 

O.lOE+01 

O.lOE-t-03 

FREQ 

339.1 

339.1 

339.3 

373.0 

379.0 

11260 

U 

O.lOE-t-Ol 

O.lOE+01 

O.lOE-fOl 

O.lOE-t-Ol 

0.10E-«-01 

O.IOE+Ol 

V 

■0.35E-06 

-0.35E-06 

-0.35E-06 

-0-.10E+08 

~0.16E+10 

O.OOE-t-lO 

W 

0.39E-08 

0.39E-06 

0.39E-04 

0.e7E>02 

0.91E-fOO 

O.llE-t-03 

0 

-O.llE-14 

-0.90E-13 

-0.92E-11 

-O.47£-f05 

-0.57E+09 

-O.34E+09 

LL 

O.lOE-03 

LK 

O.lOE-07 

O.lOE-05 

O.lOE-03 

O.lOE-01 

O.IOE+Ol 

O.lOE+03 

FREQ 

339.1 

339.1 

339.3 

373.0 

379.0 

11260 

U 

O.lOE+01 

O.lOE+01 

O.lOE+01 

O.lOE+01 

O.IOE+Ol 

O.lOE+01 

V 

-0.31E-04 

-0.31E-04 

-0.31E-04 

-O.lOE-t-06 

-0.16E+08 

0.81E4-08 

W 

0.39E-08 

0.39E-06 

0.39E-04 

0.87E-02 

0.91E+00 

O.llE+03 

0 

-0.74E-13 

-0.74E-11 

-0.75E-09 

-O.47E+03 

-0.58E+07 

-O.34E+07 

LL 

O.lOE-01 

LK 

0,10E-07 

O.lOE-05 

O.lOE-03 

O.lOE-01 

O.lOE+01 

O.lOE+03 

FREQ 

339.1 

339.1 

339.3 

373.0 

379.0 

11260 

U 

O.lOE+01 

O.lOE+01 

O.lOE+01 

O.lOE+01 

O.lOE+01 

O.IOE+Ol 

V  -0.31E-02 

-0.31E-02 

-0.31E-02 

-O.lOE+04 

-0.16E+06 

0.81E+06 

w 

0.39E-08 

0.39E-06 

0.39E-04 

0.87E-02 

0.91E+00 

O.llE+03 

0 

-0.74E-11 

-0.74E-09 

-0.75E-07 

-O.47E+01 

-0.58E+05 

-O.34E+05 

LL 

O.lOE+01 

LK 

O.lOE-07 

O.lOE-05 

O.lOE-03 

O.lOE-01 

O.lOE+01 

O.lOE+03 

FREQ 

366.9 

366.9 

367.1 

393.1 

436.6 

11310 

U 

O.lOE+01 

O.lOE+01 

O.lOE-t-01 

O.lOE-t-01 

O.lOE+01 

O.lOE+01 

V 

-0.38E+00 

-0.38E+00 

-0.39E+00 

-O.llE+02 

-0.16E+04 

0.81E-f04 

w 

0.42E-08 

0.42E-06 

0.42E-04 

0.87E-02 

0.91E+00 

0.11E+03 

0 

-0.14E-08 

-0.14E-06 

-0.14E-04 

-O.56E-01 

-0.54E+03 

-O.34E+03 

LL 

O.lOE+03 

LK 

O.lOE-07 

O.lOE-05 

O.lOE-03 

O.lOE-01 

O.lOE+01 

O.lOE-i-03 

FREQ 

16160 

16160 

16160 

16160 

16790 

38530 

U 

0.10E4-01 

O.lOE+01 

O.lOE+01 

O.lOE-fOl 

O.lOE-fOl 

O.lOE-fOl 

V 

0.44E-»'03 

0.44E-«-03 

-0.27E-01 

0.50E403 

-0.79E+02 

0.69E<»-02 

W 

-0.23E-07 

-0.23E-05 

-0.20E-03 

>0.26E-01 

0.92E4'00 

0.94E+02 

0 

0.29E'05 

0.29E-03 

0.23E-04 

O.33E<*-01 

-0.52E<»^02 

0.92E-I-01 

Tabic  7*4:  EfTects  Non-^rcular  Crocs  Section  On  the  Fundamental 
Frequency  and  Mode  Shape  (a/b=l.l,  ^=90’ ) 
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CHAPTER  8 


CONCLUSIONS 

The  primary  purpose  of  the  work  presented  in  this  report  was  to  examine  the 
effects  of  curvature  in  three  dimensional  space  and  of  out-of-roundness  on  the  vibrations  of 
rods  with  special  attention  to  rods  which  are  nearly  straight. 

Although  the  derived  governing  equations  apply  to  a  rod  with  general  space 
curvature,  the  mode  shapes  and  natural  frequencies  of  a  free  vibration  curve  were 
calculated  for  a  rod  whose  centroidal  axis  was  a  helix.  For  this  case,  both  the  curvature  and 
torsion  were  constants  rather  than  position  dependent.  Natural  frequencies  and  mode 
shapes  were  obtained  using  three  term  series  with  the  Raleigh-Ritz  method. 

It  was  found  that  the  introduction  of  even  slight  space  curvature,  as  characterized 
by  the  curvature  and  torsion  parameters,  introduced  significant  changes  on  rod  natural 
frequencies  and  especially  on  mode  shapes  compared  to  those  of  a  straight  rod.  Increasing 
either  curvature  parameter  always  increased  fundamental  natural  frequency  zind  degree  of 
coupling  among  the  displacement  variables  in  the  natural  modes. 

The  effect  of  out-of-roundness  was  examined  by  calculating  the  mode  shapes  and 
natural  frequencies  of  a  rod  of  ellipitical  cross  sections  and  varying  the  orientation  of  the 
principal  axes  of  the  ellipse  with  respect  to  the  natural  coordinates  of  space  curvature  of  the 
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rod  center  line.  It  was  found  that  the  orientation  of  the  cross  section  had  a  pronunced  effect 
on  the  coupling  of  the  displacement  variables  on  any  given  natural  mode  of  free  vibration, 
and  that  out-of-roundness  changed  the  manner  in  which  changes  of  the  torsion  and 
curvature  parameters  affected  the  rod's  natural  frequencies  and  mode  shapes. 

The  mode  shape  and  natural  frequencies  were  calculated  by  the  Raleigh-Ritz 
methods  which,  in  turn,  required  use  of  the  strain  and  kinetic  energy  functions  for  a  helical 
rod.  These  were  obtained  as  a  special  case  for  the  strain  energy  and  kinetic  energy  of  a  rod 
with  general  space  curvature  and  an  arbitrary  cross  section  shape  which  was  first  derived. 
The  general  strain  and  kinetic  energy  expressions  were  aJso  used  to  formulate  the  governing 
equations  and  natural  boundciry  conditions  of  a  rod  with  general  space  curvature.  These 
equations  were  found  to  differ  in  certain  respects  from  similiar  equations  obtained  earlier  by 
Kingsbury.  It  was  found  that  the  latter  equations  did  not  yield  a  symmetric  stiffness 
matrix  when  a  natural  mode  solution  was  introduced  but  the  newly  derived  ones  did. 

The  derived  equations  were  also  compared  to  those  obtained  by  Love  10  for  a 
circular  ring.  Imposing  the  condition  of  inextensibility  of  the  centroidal  axis  for  a  plane  ring 
yielded  the  same  in-plane  vibration  equation  as  given  by  Love.  In  the  out-of-plane 
vibration  equation  set.  Love's  equations  of  motion  were  identical  with  Kingsbury’s,  but 
slightly  different  from  those  derived  in  this  report.  As  a  further  check  on  the  derived 
energy  expression,  a  Raleigh-Ritz  solution  was  formulated  for  a  partial  circular  ring  and 
compared  with  results  presented  by  Den  Hartog  who  employed  Love’s  energy  function.  It 
was  found  that  the  energy  expression  derived  in  this  report  yields  frequencies  slightly  higher 
than  Den  Hartog’s,  especially  as  the  arc  angle  of  the  ring  became  larger.  There  is  no  way. 
however,  of  assessing  the  relative  accuracy  of  the  different  energy  expressions  from  the 
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APPENDIX  A 


STRAIN  ENERGY  FUNCTION 

Pi.  Pj.  P3  and  P4  appearing  in  equation  (2.13)  are  shown  below. 

Pi  =  Ekv"w'  ~  EkXu'w'  —  EkX'uw'  -  Ek^ow'  —  Bk  'v"w 

—  —  Ek  ' \u  w  —  Ek  'uvj  —  Gk^Xuvd  —  Gk^o^w 

—  E^:K'ou•-  EXv'v"  -  ex' —  Eu  ' v"  ~  Gk^Xw' 

—  EX^u'v'  —  Gk^u'v'  -  EXX'uv'  —  GkXov'  —  EkXov 
~  EXX  V.  V ~  EX  '^uv~  Gk^X^‘uv~  GkXo'v  —  EkX'ov 

—  EXu'u"  —  ex' um"  -  Ekou"  —  Gk^Xuv,'  —  Gko' u 

—  GkX'ou 

P2  =  ^  :  Ev""^  ~  Gk\'^  ~  EX\'^  -  EX  ~  Gk^xK^ 
-Go'^-  GA^o^ 

—  £Au  D  —  EX  xiv"  —  Ekov"  —  Gk^Xhv'  —  Gko' v 

—  EXX  mu'  —  EkXou'  —  GKXdt'u  —  BkX'^u 
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P,  =  -  -  Ek  -  Gk*w^  -  EX^v'^  -  Ex'^v^ 

2 

-  GK'^Xh'^  -  Eu"^  -  Gk\'^  -  Go '2  -  GA^o' ; 

—  Ekk'ww'  -  EkXv'w'  —  EkX'vw'  —  Eku'w'  —  Ek^Xv 

-  En^x'vw-  Gk^Xwj—  Ek.^u"w->-  Gk^uw—  Gk^X<Pw 
—  EXX^vv'  -  EXu"v'~  EX'u"v-  Gk^Xm'v-  GkX^4>v 
—  GkXou' 

P,  Ew'^  —  Ek^u^  -  Ekuw' 

2 

where 


n; 


APPENDIX  B 


MATRICES  MA  AND  MB 

In  Appendices  B.  C.  and  D.  the  following  dimensionless  parameters  have  been 
introduced. 


pi  =  r 


lk  =  /« 


11  =  a 
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B.l  Matrix  Ma 

=  Mo(y.«)  for  i=3  (ij  =  l . 12) 

mad  ,  1 )  *  (16><co28*pi'»*4+4<‘co20*ge*lk**2*pi**2)  *xy+4*co20*ll**2*pi 
1  =i‘^2+(co28+l)*ge‘*‘lk“*2=»ll**2+{co28+l)^<la*lk**2 

mad, 2)  =  ge*lk**2>^ll**2+la*lk**2 

mad, 3)  =  ge*lk**2*ll**2+la»lk*»2 

mad  ,4)  =  d6‘*co28*pi**4+(4*co20*ge*lk»*2-4*co20*ll**2)*pi**2+(-co 
1  28-l)>"ge>'lk**2*ll^>"2)-“X2 

maCl.B)  =  -ge“lk'<*2»'ll*>^2’»xz 

mad, 6)  =  -g€xlk*=^2“ll’^''2>‘x2 

mad, 7)  =  (2*col2*lk*ll*pi'‘*2-co4^ge»lk**3*ll)*X2 

rnaCl  ,8)  =  (8*co28*lk‘'pi“*3+2*co20*ge*lk*’*'3*pi)*xy+2»co28‘>la*lk>^pi 

mad, 9)  =  (6*co23*lk"ll‘"pi>^»2-co29-ge*lk*«3xll)»X2 

rnaCl,  10)  =  ((2>‘col2*ge-4>'co4)>‘lk>^pi*»2-co4*ge*lk‘‘ll**2)*xz 

rnaCl,!!)  =  2'^co20*ge*lk*ll>‘pi*xy+(-2>‘Co28*ge-2*co20)*lk*ll*pi 

mad,12)  =  ((6>‘Co23»g€-4*co29)*lk*pi**2-co29*ge*lk=^ll**2)*xz 

ma(2 ,2)  =  (256^co50*pi**4+16«co46*ge>‘lk»*2*pi**2)*xy+16'i'co46*ll*>^2 
1  >'pix>'2+(co50+l)  *‘ge>‘lk’‘*'2’^ll**2+(co50+l)*la*lk**2 

ma(2,3)  =  ge*lkxx2xllxx2+la*lkx*2 

ma(2,4)  =  -gexlk**2*ll**2*x2 

ma(2,5)  =  (256*co50*pi**4+d6*co46*ge*lk**2-16*co46*ll**2)*pi**2+( 
1  -co50-l)*ge*lk**2*ll**2)*xz 

ma(2,6)  =  -ge*lk»*2*ll**2*xz 

ma(2,7)  =  (4*col6*lk*ll*pi**2-co8*g«*lk*^3*ll)*x2 
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ma(2,8)  =  0 

ma(2,9)  =  (12‘<co-»2>‘lk“'ll*pi*'~2-co38'‘ge*lk*»3'*ll)>'xz 
ma(2,10)  =  ( (4*col6'i'ge-16‘‘Co8)>^lk*pi**2-co8*‘ge*lk*ll**2)*xz 
ma(2,ll)  *  0 

ma(2,12)  =  ( (12*co42*ge-16*co38)*lk*pi**2-co38*ge*lk*ll**2)*xz 

!na(3,3)  =  (1 296*0  o55»‘pi**4+36*co53*ge»lk**2*pi**2)*xy+36*co53*ll^* 
1  2*pi**2+(co55+l)*ge*lk»‘^'2*ll**2+(co55+l)»la*lk*»2 

ma(3,4)  =  -ge*lk**2*ll**2*xz 

ma(3,5)  =  -g€*lk**2*ll**2*x2 

ma(3.6)  =  (1296*co55*pi**4+(36*co53*ge‘‘lk**2  .  ;o53“ll**2)  *pi**2-'’ 

1  (-co55-l)*ge*lk**2*ll**2)*X2 

ma(3.7)  =  (6*col8>i‘lk*ll*pi**2-col0*ge*lk*>‘3*ll)*xz 

ma(3,8)  =  0 

ma(3,9)  -  (18*co44*lk“ll*pi*>'2-co40*ge»lk*>»3*ll)*xz 
ma(3,10)  =  ((6*col8*ge-36*col0)*lk*pi**2-col0*ge*lk*ll**2)“xz 
ma ( 3 , 1 1 )  =  0 

ma(3,12)  =  ( (18*co44*ge-36*co40)*lk*pi**2-co40*ge*lk*ll>‘*2)*xz 

ma(4 ,4)  =  (4*co20*ll*«2*pi**2+(co28+l)*ge*lk**2*ll**2)*xy+16*co28* 
1  pi**4+4*co20*ge“lk**2*pi**2 

ma(4.5)  =  ge*lk**2*ll**2*xy 

ina(4,6)  =  ge*lk**2*ll**2’!'xy 

ma(4.7)  =  (co4*ge*lk**3*ll-2*col2*lk*ll*pi**2)’*xy 
ina(4,8)  =  (8*co28*lk*pi»*3+2*co20*ge*lk**3*pi)*xz 
ma(4.9)  =  (co29*g€*lk**3*ll-6*co23*lk*ll*pi**2)*xy 
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maC^i.lO)  =  co*i>'ge“lk^ll~*2"X)'+(2>'Col2*ge-4xco4)*lk>*pi*>^2 

nia(4,ll)  =  ((2>'co28+2*co20)*g6+2*co20)>’lk*ll*pi>*X2 

ni&(‘i,12)  =  co29*ge»'lk*ll**2*xy+(6*co23*ge-4>*co29)*lk*pi*>^2 

ma(5,5)  =  (16»'co46*ll**2*‘pi**2+(co50+l)*ge*lk**2'i‘ll’i'*2)*xy+256»co5 
1  0*pi*‘^4+16*co46*ge*lk**'2*pi**2 

ma(5,6)  =  ge*lk>**2>^ll**2*xy 

ma(5.7)  =  (co8’^ge*lk“'*3»ll-4*col6>'lk*ll*pi**2)*xy 
ma(5,8)  =  0 

ma(5,9)  =  (co38>'ge^lk>'’‘3>‘ll-12xco42*lk*ll^pi**2)*xy 
ma(5,10)  =  co8xge‘‘lk»ll**2>'xy+(4*‘col6*ge-16*co8)‘*lk*pi'*'>*2 
ma(5,ll)  =  0 

ma(5,12)  =  co38>‘ge»lk“ll»-2*xy+(l2=^co42*ge-16*co38)*lk‘i‘pi“-2 

ma(6 ,6)  =  ( 36^co53'' 11'^  *2>‘pi>‘^2+(co55+l)*ge>‘lk>‘»'2>‘ 11**2)  “xy+1296*co 
1  55>'pi''*4*36«co53''ge*lk*>'2>'pi**2 

ma(6,7)  =  (colO»g«»lk**3*ll-6*col8*lk*ll>'pi**2)  *xy 

inai;6,8)  =  0 

ma(6,9)  =  (co40>‘ge“lk**3'<ll-18*co44*lk*ll*pi**2)*xy 
ma(6,10)  =  col0>‘gexlk*ll**2*xy*(6*col8>‘g€-36*col0)*lk*pi*“'2 
ma(6,ll)  =  0 

ma(6,12)  =  co40*ge*lk*ll**2*xy+(18*co44*ge-36*co40)*lk*pi**2 
ma(7,7)  =  (coll*lk**2*pi**2+col*ge*lk**4)*xy+coll*la*pi**2 
ma(7.8)  =  0 
ma(7,9)  =  0 

ma(7,10)  =  col*ge*lk**2*ll*xy 
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ma(7,ll)  =  (2*co4>'ge*col2) 
ma(7.12)  =  0 

ma(8,8)  =  (•4>‘co28“lkx=^2*pi*>«2+co20*ge*lk'‘^4)  *xy+4*co28*la*pi‘<*2 
ma(8.9)  =  0 

ma(8.10)  =  (col2>‘ge-2>'Co4)  *lk**2’‘pi»xz 
ma(8.11)  =  co20^ge*lk*’'2*ll“^xy 
ma(8,12)  =  (3“Co23*ge-2'^co29)*lk**2»pi*xz 

ma(9,9)  =  (9>'Co41>‘lkx*2>'pi**2+co3B>‘ge*lk>'»4)  *xy+9*co41*la*pi**2 
ma(9,10)  =  0 

jna(9,ll)  =  (2xco29xge-3xco23)*lkx*2xpi*xz 
ma(9,12)  =  co35xgexlk*x2xllxxy 

ma(10 , 10)  =  -lm<iax(-colxx}'-col)  +  (collxgexpi'*'»2*colxlmda+col>‘gexllx 

1  x2)  xxj'+collxgexpi»*2-colxlmda*col*^ge*llxx2+col‘‘lkx*2 

madO.ll)  =  0 

ma(10,12)  =  0 

ma(ll .  11)  =  -Im<iax(-co20xxy-co20)'*-(4*co28*ge*pi»‘*2-co20*lmda+co20x 

1  ge“llx  *2)  *xy*4xco28xgexpi**<2-co20’^lmda+co20’'‘ge*ll**2+co20*lk*x2 

ma(ll,12)  =  0 


ma(12 , 12)  =  -Imda*  (-co35“'xy-co35)  +  (9*co41*ge*pi**2-co35*iinda+co35* 
1  ge*ll*  >^2)  xxy+9xco41*g€xpi**2-co35*lmda+co35>i‘ge*ll**2+co35*lk**2 
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B.2  Matrix  Mb 


=  MfcO’.i)  fo-  1=5'  (tj  =  l . 12) 

mbd.l)  =  -  (-C028-1) 

mb (1,2)  =  la 

mb(l ,3)  =  la 

mb(l ,4)  =  0 

mb(l ,5)  =  0 

mb(l,6)  =  0 

mb(l,7)  =  0 

mb(1.8)  =  0 

mb(l ,9)  =  0 

mb(l . 10)  =  0 

mb  ( 1 , 1 1 )  =  0 

mb(l.l2)  =  0 

mb(2,2)  =  -  (-coBO-d^la 

mb(2,3)  =  la 

mb(2.4)  =  0 

mb(2.5)  *  0 

mb(2.6)  =  0 

mb(2.7)  =  0 

mb(2.8)  =  0 


mb(2.9)  =  0 
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mb(2.10)  =  0 

mb(2.11)  =  0 

mb(2.12)  =  0 

mb(3.3)  =  -(-co55-l)“la 

mb(3.‘i)  =  0 

mb(3.5)  =  0 

mb(3.6)  =  0 

mb (3. 7)  =  0 

mb(3.8)  =  0 

mb(3.9)  =  0 

mb(3.10)  =  0 

mb(3.11)  =  0 

mb(3.l2)  =  0 

mb(-;,-i)  =  -  (-C028-1) 'la 

mb(-i.5)  =  la 

mb(-i,6)  =  la 

mb  (-1,7)  =  0 

mb(4.8)  =  0 

mb (4, 9)  =  0 

mb(4,10)  =  0 

mb(4.11)  =  0 

mb(4,12)  =  0 


mb(5,5)  =  -(-coSO-D-la 
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mb(5.6)  =  la 

mb(5.7)  =  0 

mb(5.8)  =  0 

mb(5.9)  =  0 

mb(5.10)  =  0 

mb(5.11)  =  0 

mb(5.12)  =  0 

mb(6.6)  =  -(-co55-l)-la 

mb(6.7)  =  0 

mb(6,8)  =  0 

mb(6,9)  =  0 

mb(6.10)  =  0 

mb(6.11)  =  0 

mb(6.12)  =  0 

mb(7,7)  =  col*la 

mb(7.8)  =  0 

mb(7.9)  =  0 

mb(7.10)  =  0 

mb(7,ll)  =  0 

mb(7.12)  =  0 

mb(8,8)  =  co20>-la 

mb(8.9)  =  0 


mb(8.10)  =  0 
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mb(8.11)  =  0 
mb(8,12)  =  0 
mb(9.9)  =  co35'la 
mb(9.10)  =  0 
mb(9.11)  =  0 
mb(9.12)  =  0 
mbdO.lO)  =  col^xy+col 

mb(lO.n)  =  0 

mb(10.12)  =  0 
mb  (11, 11)  =  co20-xi'+co20 
mb(11.12)  =  0 
mb(12,12)  = 

B.3  Integrated  Coefficients 

The  integrated  coefficients,  col  to  co55.  shown  in  matrices  and  .\Ij  are  the 
results  of  the  product  of  functions  ^  and  rj  integrated  from  =  0  to  ^ 

Functions  ^  a:  d  rj  as  well  as  the  relations  between  con  (n  =  1 . 55)  and  functions 

.  T]  are  given  in  table  B-1. 

From  table  B-1,  we  arrive  at  the  following  values. 

col=  0.5 
co2=  0 
co3=  0 
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co;=  -  0.2122065907891938 
co5=  0 
co6-  0 
co7=  0 

co8=  -  0.042-i4131815783876 
co9»  0 

col0=-  0.01818913635335947 
coll=  O.b 

col2=  0.4244131815783876 
col3=  0 
col4=  0 
col5=  0 

col6=  0.169765272631355 
col7=  0 

col8=  0.1091348181201568 

col9=  0 

co20=  0.5 

co21=  0 

co22=  0 

co23=-  0.2546479089470326 

co24=  0 

co25=  0 

co26=  0 

co27=  0 

co28=  0.5 

co29=  0.3819718634205488 

co30=  0 

co31*  0 

co32=  0 

co33=  0 

co34=  0 

co35=  0.5 

co36=  0 

co37=  0 

co38=-  0.272837045300392 
co39=  0 

c o40= -  0 . 0707355302630646 
co41=  0.5 

co42=  0.3637827270671894 

co43=  0 

co44=  0.1414710605261292 

co45=  0 

co46=  0.5 

co47=  0 

co48=  0 


co49=  0 

co50=  0.5 

co5l=  0 
co52=  0 

co53=  0.5 

co54=  0 

cc55=  0.5 


APPENDIX  C 


MATRICES  MC  AND  MD 


C.l  Matrix  Me 


M,  (!j')  =  M,  (y.i)  for  {ij=l . 8) 

mc(l  ,1)  =  ( ( (2*co70-2>'co65-2»co61-*-2*co59)*ge*kl»kkl»*2'*-(-2*co73-2» 

1  co69+-.'Co66-2“Co64+2-co58)‘‘ge»kl'‘kkl  +  (2»co70-2*co65+2*co61-2*co 

2  59)  “ge-kl)  xlk-*2-*-  (2x0  070- 2*co65+2*co61-2*co59)  *kl-^*3*kkl**2+(-2 

3  X<;o73-2xco69+2*co64-4xco60-2xco58)*kl*»3*kkl+(2*co70+2*co65+2*c 

4  o61+2xco59) xklxx3) xiixx2+ ( ( (co73-2“Co66+co64)*ge*kl**2*kkl**2+( 

5  - 2^0 070+2*0 o65-2xco61+2xco59)*gexkl**2*kkl+(co69+2*co60+co58)xg 

6  e*kl* -2) X  Ik* *2+ (0  069+2*0060+0 o58)*kl*x4+kkl **2+( *2*0070-2*0065- 

7  2xoo61-2xoo59)xklxx4xkkl+(oo73+2xoo66+oo64)*kl**4)*xy+( ( (oo69-2 

8  xoo60+oo58)xgexkklxx2+(-2xoo70+2»oo65+2*oo61-2*oo59)*ge*kkl+(oo 

9  73-2xoo66+oo64)“ge) xlkxx2+(oo73-2*co66+co64)*kl**2*kkl**2+(-2*c 
o70+2>'co65-2xoo6l+2*oo59)xkl*x2*kkl+(oo69+2*oo60+oo58)*kl**2) *1 

;  1*  2+ ( (o 069-2*0060+0058) *kkl**2+(-2*co70+2*oo65+2*oo61-2*oo59) * 

<  kkl+oo73-2*oo66+co64) *la*lk**2 

mo (1,2)  =  (((( (oo54-oo52-oo39+co37)*ge*k2+(oo48-oo46-oo31+oo29}*g€ 

1  *kl) *kkl+ (-oo55+co53+oo40-oo38)*ge*k2+(-oo47-co45+oo30+oo28)*ge 

2  *kl) *kk2+ ( (-oo47+oo45-oo30+oo28)*ge*k2+(-oo55+oo53+oo40-oo38) *g 

3  e*kl) *kkl+ (oo48-co46+co31-oo29)*ge*k2+(co54+oo52-oo39-co37) *ge* 

4  kl)*lk**2+( ( (oo48-oo46+oo31-co29)*kl*k2**2+(oo64+oo52-co39-co37 

5  ) *kl**2»k2)*kkl+(-co47-oo45-oo30-co28)*kl*k2**2+(-oo55'Oo53+oo4 

6  0+oo38)*kl**2*k2) *kk2+( (-oo55+oo53-oo40+co38)*kl*k2**2+(-oo47-c 

7  o45-oo30-co28)*kl**2*k2)*kkl+(oo54+oo52+oo39+oo37)*kl*k2**2+(oo 

8  48+oo46+oo31+oo29)*kl**2*k2)*ll*X2+{( ( (oo55-oo53-oo40+oo38) *ge* 

9  kl*k2*kkl+(-oo54-oo52+oo39+oo37)*ge*kl*k2)*kk2+(-oo48+oo46-oo31 
+oo29) *ge*kl*k2*kkl+ (oo47+oo45+oo30+oo28) *g«*kl*k2) *lk**2+ ( (oo4 

;  7+00'  5+oo30+oo28)*kl**2*k2**2*kkl  +  (-co48-co46-co31-oo29)*kl**2* 

<  k2**2)*kk2+(-oo54-co52-oo39-oo37)*kl**2*k2**2*kkl+(oo55+oo53+oo 

=  40+oo38)*kl**2*k2**2)*xy+((((oo47-oo45-oo30+co28)*ge*kkl+(-oo48 

>  +oo46+oo31-oo29)*ge)*kk2+(-oo54+co52+oo39-co37)*ge*kkl+(oo55-oo 

?  53-oo40+oo38) *ge) *lk**2+( (oo55-co53-oo40+oo38)*kl*k2*kkl+(-oo54 
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-co52+co39+co37)  *kl*k2)  >'kk2+(-co48+co46-co31+co29)  >‘klxk2*kkl  +  (c 

1  o47+co45+co30+co28) *kl>'k2)*ll**2+( ( (co47-co45-co30+co28) *kkl-co 

2  48+co46+co31 -co29)*kk2+(-co54+co52+co39-co37) ^kkl+co55-co53-co4 

3  0+co38)^la-lk^'“<2 

mcd  ,3)  =  ((((-co69+2*co60-co58)*ge«kkl**2+{2*co70-2*co65-2*co61+2 

1  *co59)*ge^kkl  +  (-co73-*-2xco66-co64)*ge)*lk**2+(-co73+2*co66-co64) 

2  *kl**2*kkl»*2+(2*co70-2*co65+2*co61-2*co59)*kl**2*kkl+(-co69-2* 

3  co60-co58)  *kl'^*2)’^ll’*^*2+(  (co73-2*co66+co64)*ge*kl**2*kkl**2+(-2 

4  *co70+2*co65-2>i'co61+2’i‘co59)  *ge»kl**2»kkl+ (€069+2*0060+0058)  *ge* 

5  kl**2)*lk*>'2+(oo69+2*oo60+oo58)*kl**4*kkl**2+( *2*0070-2+0065-2* 

6  oo61-2»oo59)*kl**4*kkl+(co73+2*co66+co64)*kl**4)*xz+(( (~oo70+oo 

7  65+oo61-oo59)*ge*kl*kkl**2+(oo73+co69-2*oo66+oo64-oo58)*ge*kl»k 

8  kl+ ( -oo70+oo65-oo61+oo59)*ge*kl)*lk**2+(-oo70+oo65-oo61+oo59) *k 

9  l**3*kkl**2+(oo73+oo69-oo64+2*oo60+oo58)*kl**3*kkl+(-oo70-oo65- 
oo61-co59) *kl^*3)*ll*xy+( ( (oo70-oo65-oo61+oo59)*ge*kl*kkl**2+ (- 

:  Co73-oo69+2*oo66-oo64+oo58) *ge*kl*kkl+(oo70-oo65+oo61-oo59) *ge* 

<  kl) “lk*»2+(co70-oo65+oo61-oo59)*kl**3*kkl**2+(-oo73-oo69+oo64-2 
=  *oo60-oo58) *klx*3*kkl+ (co70+oo65+oo61+oo59)*kl**3)*ll 


mo  (1 .4)  =  ( ( ( ( (-oo47+oo45+oo30-oo28)*ge*kkl+(oo48-oo46-oo31+oo29) * 

1  ge) xkk2+(oo54-oo52-oo39+oo37)*ge*kkl+(-oo55+oo53+oo40-oo38)>'ge) 

2  *lkx*2+( (-oo55+oo53+oo40-oo38)xkl*k2*kkl+(oo54+oo52-oo39-oo37)* 

3  kl*k2)xkk2*(oo48-oo46+oo31-oo29)xkl*k2*kkl+( -0047-0 o45-oo30-oo2 

4  8)-klxk2)*llxx2+(((oo55-oo53-oo40+oo38)*ge*kl*k2*kkl+(-oo54-oo5 

5  2+co39+oo37)xgexkl*k2)*kk2+(-oo48+oo46-oo3l+oo29)*ge*kl*k2»kkl+ 

6  (oo47+co45+oo30+oo28) *ge*kl*k2)»lk**2+((oo47+oo45+oo30+co28)*kl 

7  **2xk2*x2*kkl+(-oo48-oo46-oo31-oo29)*kl**2*k2**2)*kk2+( -0054-00 

8  52-co39-co37)*kl**2*k2x*2*kkl+(oo55+oo53+oo40+oo38)*kl**2*k2**2 

9  )*X2+( ( ( (-Co48+oo46+oo3l-co29)xge*kl*kkl+(oo47+oo45-oo30-oo28)* 
ge^kl) xkk2+ (oo55-oo53-co40+oo38)“ge*kl*kkl+ (-oo54-oo52+co39+oo3 

;  7)xge*kl) *lk**2+( (-oo54-oo52+oo39+oo37)*kl**2*k2*kkl+(oo55+oo53 
<  -oo40-oo38) *kl**2*k2) *kk2+ (oo47+co45+oo30+oo28)*kl**2*k2*kkl+(- 

=  co48-co46-co31 -oo29)*kl**2*k2)*ll*xy+( { ( (co54-oo52-oo39+oo37) *g 

>  e*k2*kkl+(-oo55+oo53+oo40-oo38)*ge*k2)*kk2+(-oo47+oo45-oo30+oo2 

?  8)*ge*k2*kkl+(oo48-oo46+oo31-co29)*g«*k2)*lk**2+((oo48-oo46+oo3 

l-oo29)*kl*k2**2*kkl+(-oo47-oo45-oo30-oo28)*kl*k2**2)*kk2+(-oo5 

1  5+oo53-oo40+co38)*kl*k2**2*kkl+(oo54+oo52+oo39+oo37)*kl*k2**2)* 

2  11 

mod ,5)  =  (((oo68-oo75)*kl*kkl+(oo72+oo63)*kl)*lk*ll*pi+((co62-co7 

1  l)*ge*kkl+(oo74-oo67)*ge)*lk**3*ll)*xz+(((-co72-co63)*kl**2*kkl 

2  +(oo75+oo68)*kl**2)*lk*pi+((oo67-co74)*ge*kl*kkl+(co71+oo62)*ge 

3  *kl)»lk**3)*xy+((oo72-oo63)*kkl-oo75+oo68)*la*lk*pi 
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me  (1 ,6)  =  ( ( (2>'col9-2*co21)’‘kl*kkl*(2“co20+2’*col8)*kl )  >‘lk*ll'^pi+(  ( 

1  co7-co9)  “ge-kkl*  (colO-co8)*ge)»lk*>‘3‘^ll)*xz+(  ( (-2*co20-2>'col8)* 

2  klx“'2>‘kkl+(2'<co2l+2“col9)=^kl**<2)*lk^pi+((co8-col0)*ge*kl*kkl+(c 

3  o9+co7)xge^kl)  xlk**3)  *xy+(  (2*co20*2*col8)’»kkl-2*co2l+2'<col9)*la 

4  *lk«pi 

mc(l  ,7)  =  (((co68-co75)*ge*kl*kkl+(co72+co63)>S‘ge*kl)*lk*pi+((co62- 

1  co71)*ge*kkl+(co74-co67)»ge)*lk<'ll**2v({co71+co62)*kl**2*kkl+(- 

2  co74-co67)‘^kl*^2)*lk)^x2+((co67-co74)*ge*kl*kkl+(co71+co62)*ge* 

3  kD’i'lk^ll'^xy+C  (co63-co72)*ge»kkl+(co75-co68)*ge)*lk*ll*pi+((co7 

4  4-co67)“kl>*kkl+(*co71-co62)*kl)*lk'^ll 

me (1 ,8)  =  ( ( (2*col9-2*co21)*ge*kl*kkl+(2*eo20+2*col8)*ge*kl)*lk*pi 

1  +( (co7-co9) *gexkkl+(colO-eo8)*ge)^lk*ll*‘*2+( (eo9+eo7) xklx*2*kkl 

2  +(-co8-colO)xklxx2) Xik)xx2+((co8-col0)xgexkl*kkl+(co9+co7)xgexk 

3  D'^lkxll^'xy+C  (2*eol8-2xeo20)xgexkkl+(2xco21-2*col9)xge)*lkxllxp 

4  i  +  (  (c  JlO-co8)xklxkkl-*'(-eo9-co7)*kl)xlk*ll 

me (2.2)  =  ( ( (2>‘eo44-2xco35-2*eo27+2xco25) *ge*k2xkk2x*2+(-2xeo5l-2x 

1  eo43+4*co36-2xco34+2«eo24)xge*k2xkk2+(2xeo44-2»eo35+2xco27-2*co 

2  25)xgexk2)xlkxx2+(2xeo44-2*eo35+2xco27-2xeo25)xk2xx3*kk2xx2+(-2 

3  Xco5l-2xeo43+2*co34-4xco26-2xeo24)*k2*x3*kk2+(2*eo44+2*eo35+2xe 

4  o27+2xco25)xk2x»3)xllxx2-''(((co51*2xeo36+co34)*gexk2xx2xkk2xx2+( 

5  -2>'eo44*2xeo35-2xeo2742*co25)xge«k2*x2xkk2+(eo43+2xco26+co24)xg 

6  e*k2xx2)xlkxx2+(eo43+2xeo26+eo24)xk2xx4*kk2*x2*(-2xco44-2xco35- 

7  2xeo27-2xeo25)xk2**4xkk2+(eo5l+2xeo36+eo34)xk2xx4)xxy+( ( (eo43-2 

8  xeo26+co24) xge»kk2xx2+(-2xco44+2xco35+2xco27-2xeo25)*gexkk2+(co 

9  51-2*eo36+eo34)xge)xlkx*2+(eo51-2»co36+co34)xk2*x2*kk2xx2+(-2xe 
o4-.+2xco35-2>'eo27+2xeo25)xk2x*2*kk2+(co43+2*eo26+co24)xk2xx2)xl 

;  lxx2*( (co43-2»co26+eo24)*kk2*x2+(-2*eo44+2»eo35+2“eo27-2»eo25)x 

<  kk2+co51-2’'eo36+eo34)xlaxlkx*2 

me  (2.3)  =  (((( (-eo47+co45+eo30-eo28)*gexkkl  +  (co48-eo46-eo31+eo29) * 

1  ge) »kk2+ (eo54-co52-eo39+co37)*ge*kkl+(-eo55+eo53+co40-co38)xge) 

2  *lkxx2+( (-Co55xeo53+co40-eo38)*kl*k2*kkl+(eo54+eo52-co39-eo37)* 

3  klxk2)*kk2+(co48-co46+eo31-co29)*kl*k2*kkl+(-co47-co45-co30-co2 

4  8)xklxk2) *ll**2+(((co55-co53-co40+co38)*ge*kl*k2*kkl+('co54-co5 

5  2+co39+co37)*ge*kl*k2)*kk2-*(-eo48+co46-eo31+co29)*ge'*kl*k2*kkl  + 

6  (eo47+eo45+co30+co28)*ge*kl*k?)*lk**2+((co47+co45+co30+eo28) *kl 

7  x*2*k2**2*kkl+(-eo48-co46-co31-co29)*kl**2*k2**2)*kk2+(-co54-eo 

8  52-co39*co37)*kl**2»k2**2*kkl+(co55+co53+co40+co38)*kl*’'‘2*k2**2 

9  )*xz*( (( (-co54+co52+co39-co37)*ge*k2*kkl+(co55-co53-co40+co38)* 
ge*k2)*kk2+(co47-co45'*-co30-co28)*ge*k2*kkl+(*co48+co46-co31-*-co2 

;  9)*gexk2)*lk**2-*-((-co48+co46-co31+eo29)*kl*k2**2*kkl+(co47+eo45 

<  +eo30+eo28)  xklxk2**2)  *kk2'*-(co65-co53+co40*co38)*kl*k2*x2*kkl+(- 


W 
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=  co54-co52-co39-co37)  >^kl'‘k2“‘“2)*ll'»xy+(  ( ( (co48-co46-co31+co29)~g 

>  e^kl'<kkl  +  (-co47-co45+co30+co28)'‘ge»kl)*kk2+(-co55+co53+co40-co3 

?  8)xge“kl>^kkl+(co54+co52-co39-co37)*ge'i‘kl)*lk**2+((co54+co52-co3 

9-co37)  >‘klx*2>‘k2'i'kkl  +  (-co55-co53+co40+co38)  *kl**2*k2)*kk2+(-co4 

1  7-co45-co30-co28)^kl**2‘^k2»^kkl+(co48'*-co46+co31+co29)*kl**2*k2)* 

2  11 

me  (2 ,4)  =  ( ( ( ( -co43+2*co26-co24) *ge*kk2»*2+ (2*co44-2*co35-2*co27+2 

1  >i'Co25)  *'ge^kk2+(-co51+2>fCo36-co34)>*^g€)*lk**2+(-co51+2*co36-co34) 

2  *k2**2*‘kk2**2+(2>‘Co44-2*co35+2*co27-2*co25)*k2**2*kk2+(-co43-2* 

3  co26-co24)*k2>^*2)*ll**2+((co5l-2*co36+co34)*ge*k2**2*kk2**2+(-2 

4  *co44+2*co35-2*co27+2*co26)*ge*k2<^^2*kk2+(co43+2*co26+co24)*ge* 

5  k2'^»'2)*lk**2+(co43+2*co26+co24)*k2«*4f‘kk2**2+(-2*co44-2*co35-2* 

6  co27-2“co25)>i‘k2**4^kk2+(co51+2*co36+co34)*k2**4)»X2+(  ((-C044+C0 

7  35+co27-co25)  »ge>‘k2*kk2*>^2+(co51+co43-2*co36+co34-co24)*ge*k2“‘k 

8  k2+ (-co44+co35-co27+co25) *ge’^k2)*lk**2+(-co44+co35-co27+co25) *k 

9  2**3^kk2<'“"2'^(co51+co43-co34+2»co26+co24)*k2**3*kk2+(-co44-co35- 
co27-co25)  *k2x‘=3)*ll>^xy+(  ( (co44-co35-co27+co25)*ge=*k2*kk2**2+(  - 

:  co51-co43+2‘‘Co36-co34+co24)»ge*k2*kk2+(co44'Co35+co27-co25)*ge* 

<  k2)*lkx"2-t-(co44-co35+co27-co25)*k2**3*kk2**2+(-co51-co43+co34-2 

=  >'co26-co24) “k2**3“kk2+(co44+co35+co27+co25)*k2**3)*ll 


me  (2 .5)  =  ( ( (eo42-eo57)*k2*kk2+(eo60+eo33)»k2)*lk*ll>i‘pi+((eo32-eo4 

1  9)»ge>'kk2+(eo56-eo41)*ge)‘^lk>'»3*ll)»'xz+(((-eo50-eo33)=^k2>‘»2“‘kk2 

2  +(eo57+eo42)  >‘k2‘‘>‘2)»lk*pi+(  (eo41-co56)*ge*k2*kk2+(co49+eo32)*ge 

3  “^k2)  >‘lk>'>‘3)*xy+(  (eo50-eo33)*‘kk2*eo57+eo42)*la‘i‘lk*pi 

me  (2. 6)  =  ( ( (2'-eol5-2>‘eol7)  »k2*kk2+(2*eol6+2*eol4)*k2)>“lk*ll*pi+ ( ( 

1  eo3-co5)  “ge>‘kk2+  (eo6-eo4)’»ge)>‘lk**3*‘ll)  *xz+  ( ( (-2*eol6-2’'eol4)  “k 

2  2«’'2'kk2-^(2*col7+2'^eol5)*k2>»*2)*lk’“pi+((eo4-eo6)*ge*k2*kk2^(eo5 

3  +eo3)  >^ge»k2)  *lk*>‘3)  •'xy+(  (2*eol6-2*eol4)’*kk2-2*eol7+2*eol5)*la“‘l 

4  kxpi 

me  (2 ,7)  =  ( ((eo42-eo57)*ge*k2*kk2•'■(co50+eo33)*ge*k2)*lk*pi+((eo32- 

1  eo49)  >^ge»kk2+  (eo56-eo41 )  *ge)  *lk*ll’^*2+  ( (eo49+eo32)  ♦k2**2*kk2+  ( - 

2  eo56-co41)*k2*«2)i‘lk)*xz+((eo41-eo56)*g**k2*kk2-t-(eo49+eo32)*ge* 

3  k2)*lk*ll*xy+((eo33-eo50)*ge*kk2+(eo57-eo42)*ge)*lk*ll*pi+((eo5 

4  6-eo41)*k2*kk2+(-eo49-eo32)*k2)*lk*ll 

me  (2 .8)  =  ( ( (2«eol5-2*eol7)  ■i‘g«*k2*kk2+  (2*eol6+2*eol4)  *ge*k2)  *lk*pi 

1  ♦((eo3-eo5)*ge*kk2+(eo6-eo4)»ge)*lk*ll**2+((eo5+co3)*k2**2*kk2+ 

2  ( -co6-eo4) ♦k2**2) ♦Ik) ♦xz+ ( (eo4-eo6) *ge*k2*kk2+(co5+eo3) *ge*k2) ♦ 

3  lk*ll*xy+( (2*eol4-2*col6)*ge*kk2+(2^eol7-2*col5)*g«)^lk*ll*pi+( 

4  (eo6-eo4) ♦k2*kk2+ ( -eo5-co3) ♦k2) *lk*ll 
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me  (3 ,3)  =  ( ( (-2xco70+2’*co65+2"Co61-2*co59)>‘ge*kl*kkl**2+(2*co73+2* 

1  co69-4*co66+2*co64-2*co58)*ge*kl*kkl+(-2*co70+2*co65-2>*co61+2*c 

2  o59)*ge'"kl)>^lk^*2+(-2>*co70+2=*co65-2*co61+2*co59)>!‘kl>‘*3'i‘kkl**2+( 

3  2>»co73+2=^co69-2^co64+4*co60+2^co58)*kl*«3*kkl+(-2*co70-2*co65-2 

4  *co6l-2*co59)*kl*‘i‘3)*ll*X2+(((co69-2*co60+co58)*ge*kkl**2+(-2>^c 

5  o70+2*co65+2*co61 -2*0059)  *ge‘i‘kkl+(co73-2*co66+co64)*ge)i'lk**2+( 

6  co73-2*co66+co64)*kl**2*kkl**2-*'(-2*co70+2*co65-2*co61+2*co59)*k 

7  I**2*kkl+(co69+2*co60+co58)*kl**2)*ll**2*xy+((co73-2*co66+co64) 

8  *ge*kl**2*kkl**2+(-2*co70+2*co65-2*co61+2*co59)*ge*kl**2*kkl+(c 

9  o69+2*co60+co58)  *ge*kl**2)*lk**2-*-(co69+2*co60+co58)*kl**4*kkl** 
2+(-2“'Co70-2*co65-2*co61-2*co59)*kl**4*kkl+(co73+2*co66+co64)*k 

;  l>i“i<4 

me (3 ,4)  =  ( ( ( ( (-eo54+eo52+eo39-eo37)*ge*k2+(-co48+eo46+co31-co29) * 

1  ge^kl) *kkl+ (co55-eo53-eo40+eo38)*ge*k2+ (eo47+eo45-co30-eo28)*ge 

2  *kl) *kk2+ ( (eo47-eo45+eo30-eo28)*ge*k2+(co55-eo53-eo40+eo38) *ge* 

3  kl)*kkl+(-co48-*'Co46-co31-*-co29)*ge*k2+(-co54-co52+co39+co37)*ge* 

4  kl) »lk>'x2+ ( ( (-eo48+eo46-eo31+eo29)*kl*k2**2+(-co54-eo52+co39+eo 

5  37) xkl*x2xk2) *kkl+(eo47+eo45+eo30+eo28)*kl*k2**2+(co55+eo53-eo4 

6  0-eo38)xkl*x2*k2) *kk2+( (eo55-eo53+co40-eo38)*kl»k2**2+(eo47+eo4 

7  5+eo30+eo28)xkl*»2xk2)*kkl+(-eo54-eo52-eo39-eo37)*kl*k2**2+(-eo 

8  48-eo46-eo31-eo29)xklx*2xk2)xllxxz+((((co47-co45-co30+co28)*ge* 

9  kkl+ ( -co48+eo46+eo31 -eo29)xge) *kk2+( -co64+eo52+eo39-co37) *ge*kk 
l+(co56-eo53-eo40+eo38) xg€)xlkxx2+((co55-eo53-eo40+eo38)*kl*k2* 
kkl+(-eo54-eo52+eo39+eo37)»kl*k2)xkk2+(-eo48+eo46-eo31+eo29)*kl 

<  xk2*kkl+(eo47+eo45+eo30+eo28)*kl*k2)xll**2*3cy'-(  ((co55-eo53-co40 

=  +eo38)xg€xklxk2xkkl+(-eo54-eo52+co39+eo37)*ge*kl*k2)*kk2+(-eo48 

>  ♦eo46-co31+eo29) Xgexkl*k2xkkl+(eo47+eo45+eo30+eo28)*ge*kl*k2)“l 

?  k*x2+( (eo47+eo45+eo30+eo28)*kl**2xk2**2*kkl+(-eo48-eo46-eo31-eo 

29)*kl*x2*k2xx2)*kk2+(-eo54-co52-co39-eo37)*kl**2*k2**2*kkl'*-(eo 
1  55+eo53+eo40+eo38) *kl**2*k2**2 

me  (3 ,5)  =  ( ( (-eo72-co63)xkl**2*kkl-*-(eo75+eo68)*kl**2)*lk*pi+(  (eo67 

1  -eo74)xge>^kl*kkl+(eo7l+eo62)xge*kl)*lk**3) *xz+( ( (eo75-eo68) *kl* 

2  kkl+ (-eo72-eo63) *kl) *lkxll*pi+((eo71-eo62)*ge*kkl+(eo67-eo74)*g 

3  e)*lk*x3*ll)*xy 

me (3 ,6)  =  ( ({-2*eo20-2*col8)*kl**2*kkl+(2*co21*2*col9)*kl**2)*lk*p 

1  i+((eo8-eol0)*ge*kl*kkl*(eo9+co7)*ge*kl)*lk**3)i'XZ+(((2*eo21-2* 

2  col9)*kl*kkl+(-2*co20-2*eol8)*kl)*lk*ll*pi+((co9-co7)*ge*kkl+(e 

3  o8-eol0)*ge)*lk**3*ll)*xy 

me (3,7)  =  ( ((co72-co63)*ge*kkl+(co68-co75)*ge)*lk*ll*pi+( ( (co67-co 

1  74) *ge-eo74+co67) *kl *kkl+ ( (co71+co62) *g«+co71+co62) »kl) *lk*ll) * 

2  xz+((eo71-eo62)*ge*kkl+(eo67-co74)*ge)*lk*ll**2*xy*’((eo68-eo75) 
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3  *g€>'kl*kkl  +  (co72+co63)  *ge''kl)’i'lk*pi  +  (  (co7l+co62)  “‘kl*=*‘2*kkl  +  (-co 

4  74-C067) •kl**2) *lk 

me  (3.8)  =  ( ( (2»co20-2“col8)>^ge>^kkl+(2*cold-2*co21)*ge)*lk*lI^pi+(  ( 

1  (co8-colO) *ge+co8-colO) *kl*kkl+( (co9+co7)*ge+co9+co7)*kl)*lk*ll 

2  )“'XZ+(  (co9-co7)=i<ge*kkl  +  (co8-colO)*ge)*lk*ll**2’*‘xy+((2*col9-2*co 

3  21 )*ge'“kl»kkl+( 2*0020+2*0018) *ge*kl)*lk*pi+((co9+co7)*kl**2*kkl 

4  +(-oo8-oolO)*kl**2)*lk 

mo (4 ,4)  =  (( (-2*oo44+2*oo35+2*oo27-2*oo25)*ge*k2*kk2**2+(2*oo51+2* 

1  oo43-4*oo36+2*oo34-2*oo24)*g€*k2*kk2+(-2*oo44+2*oo35-2*oo27+2*o 

2  o25)*ge*k2)*lk**2+(-2*oo44+2»oo35-2*oo27+2*oo25)*k2**3*kk2**2+( 

3  2*oo51+2“oo43-2*oo34+4*co26+2*oo24) *k2**3*kk2+ ( -2*oo44-2*oo35-2 

4  *co27-2*oo25)*k2**3) *ll*xz+( ( (oo43-2*oo26+oo24)*ge*kk2**2+(-2*o 

5  o44+2*co35+2*co27-2*oo25) *ge*kk2+(oo51-2*oo36+oo34) *ge) *lk**2+ ( 

6  0 o5 I - 2*0 o36+oo34)*k2**2*kk2**2+ (-2*0044 +2*oo35-2*co27+2*co25)*k 

7  2**2*kk2+(co43+2*oo26+oo24)*k2**2)*ll*»2*xy+( (oo5l-2*oo36+oo34) 

8  *ge*k2**2*kk2**2+(-2*oo44+2*co35-2*oo27+2*oo25)*ge*k2**2*kk2+(c 

9  o43+2*oo26+co24) *ge*k2**2)*lk**2+(oo43+2*oo26+oo24)*k2**4*kk2** 
2+(-2*co44-2*co35-2*oo27-2»oo25)*k2**4*kk2+(oo51+2*oo36+co34)*k 

;  2“  *4 

mo  (4 .5)  =  (( (-oo50-oo33)*k2**2*kk2+(co57+oo42)*k2**2)*lk>‘pi+(  (oo41 

1  -co66) *ge*k2*kk2+(co49+co32)>‘ge*k2)*lk**3)*xz+( ((oo57-oo42)*k2» 

2  kk2+(-co50-oo33)*k2)>'lk*ll*pi+((oo49-oo32)*ge*kk2+(co41-oo56)»g 

3  e)*lk**3*ll)*xy 

mo (4 ,6)  =  ( ( (-2*0016-2*0014) *k2**2*kk2+ (2*0 ol7+2*col5)*k2* *2) “Ik-p 

1  i+( (oo4-oo6)*ge*k2*kk2+(oo5+oo3)*ge*k2)*lk**3)*xz+( ( (2*0017-2*0 

2  ol5)*k2*kk2+(-2*ool6-2*ool4)*k2)*lk*ll*pi+( (oo5-oo3)*ge*kk2+(oo 

3  4-oo6) *ge) *lk**3*ll) *xy 

mo (4 ,7)  =  ( ( (oo50-oo33) *ge*kk2+(oo42-oo57)*ge)*lk*ll*pi+( ( (oo41-oo 

1  56) *ge-oo56+oo41) *k2*kk2+( (oo49+co32)*ge+oo49+oo32)*k2) *lk*ll) * 

2  xz+( (oo49-oo32)*ge*kk2+(oo41-oo56)*ge)*lk*ll**2*xy+((oo42-oo57) 

3  *ge*k2*kk2+(oo50+oo33)*g«*k2)*lk*pi+((oo49+oo32)*k2**2*kk2+(-oo 

4  56-oo41)*k2**2)*lk 

mo (4,8)  =  ( ( (2*ool6-2*ool4)*ge*kk2+(2*ool5-2*ool7)*ge)*lk*ll*pi+( ( 

1  (oo4-oo6)*ge-oo6+oo4)*k2*kk2+((oo5+oo3)*ge+oo5+oo3)*k2) *lk*ll)* 

2  xz+((oo5-oo3)*ge*kk2+(oo4-oo6)*ge)  Ik*ll**2*xy+((2*ool5-2*ool7) 

3  *ge*k2*kk2+ (2*0016+2*0014) *ge*k2)*lk*pi+( (oo5+oo3) *k2**2*kk2+(- 

4  oo6-oo4)*k2**2)*lk 

mo(5,5)  =  (oo78*lk**2“pi**2+oo76*ge*lk**4)*xy+oo78*la*pi**2 


81 


me  (5,6)  =  (2~co23-lk^^2~pi“^2-*-coll>'ge‘‘lk>‘>^4)*xy+2*co23'<la*pi>“i‘2 

mc(5,7)  =  (co77*ge-co77)^lk'‘“2>'pi’'X2+co76xge*lk*>"2“ll>'3cy 

me  (5,8)  =  (2*<eo22*ge-eol2)>'lk>‘^2»'pi>‘Xz+coll*ge*lk**2*ll*xy 

me  (6,6)  =  (■i*col3*lk‘*>‘2’^pi*>^2->-col»ge’^lk*^4)*xy*4*col3’^la*pi**2 

me(6,7)  =  (eol2>^ge-2>“eo22)><lk>'“2'^pi*xz+coll*ge*lk**2*ll>!'xy 

me(6,8)  =  (2''eo2^ge-2>‘co2)"‘lk*«2*pi*xz+col*ge*lk‘f‘*2*ll*xy 

me  (7 , 7)  =  -Imda*  (-eo76*xy-eo76)  +  (co78'i‘ge*pi**2-co76*lmda+co76*ge^l 
1  l''*2)  >‘xy+co78*ge''pi*»^2-eo76>‘lmda+co76>‘ge*ll**2+co76*lk*‘*2 

me  (7 ,8)  =  -Imda*  (-eoll>^xy-coll)  +  (2*co23‘^ge*pi*>!‘2-coll>‘lmda+coll»ge 
1  2)  >‘xy*2>-eo23*ge>'pi«*2-eoll*lmda+eoll’''ge*ll**2+coll*lk*^2 

me  (8,8)  =  -lmdax(-eol*'xy-eol)  +  (4*col3>*ge*pi**2-col*lmda+col*ge*ll>‘ 

1  *2)  'col3xgexpix»2-eolxlmda+colxgexll**2+eol*lk*x2 

C.2  Matrix  Md 

Mj(ij)  =  MjO'.il  for  (ij=l . 12) 

md(l ,1)  =  - ( (-eo69+2xeo60-eo58)xkkl**2+(2xco70-2xco65-2»co61+2xeo5 

1  9)  >:kkl-co73+2xco66-co64)  “^la 

md(l ,2)  =  - ( ( (-eo47+eo45*eo30-eo28)xkkl+co48-co46-eo31+eo29) xkk2+( 
1  co54-eo52-eo39+co37) xkki-co55+eo53+co40-eo38)xla 

md(l,3)  =  0 

md ( 1 , 4 )  =0 

md(1.5)  =  0 

md(l,6)  =  0 

md(l,7)  =  0 


md(l ,8)  =  0 
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1 


1 


1 


md(2,2)  =  -((-co.;3+2*co26-co24)-kk2*-2+(2«co44-2>‘co35-2>‘Co27+2>‘co2 
5)  “kk2-co51+2''Co36-co34)  -  la 

ind(2.3)  =  0 

md(2.4)  =  0 

md(2,5)  =  0 

md(2.6)  =  0 

md(2.7)  =  0 

md(2.8J  =  0 

md(3.3)  =  - ( (-co69+2-co60-co58)-kkl*-2+(2-co70-2»co65-2*co61+2’*co5 
9) -kkl -co73+2-co66-co64) »la 

md(3 ,4)  =  - ( ( ( -co47+co45+co30-co28)*kkl+co48-co46-co31+co29) -kk2+( 
co54-co52- V  ,'39+co37)  -kkl  -co55+co53+co40-co38)  *la 


md(3.5)  =  0 
md(3.6)  =  0 
md(3.7)  =  0 
md(3.8)  =  0 

md(4 ,4)  =  - ( (-co43+2-co26-co24)*kk2**2+(2»co44-2-co35-2»co27+2-co2 
1  5) -kk2-co51+2-co36-co34) -la 


md(4,5) 

=  0 

nid(4.6) 

=  0 

md(4.7) 

=  0 

md(4.8) 

=  0 

md(5.5) 

=  co76*la 

md(5,6) 

=  coll*la 

md(5.7)  =  0 
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md(5 , 8) 
md(6,6) 
md(6 ,7) 
md(6.8) 
md(7.7) 
md(7.8) 
md(8 , 8) 


0 


col -la 


0 

0 

co76*xy»co76 
col  l*xy  •♦•coll 
col■'xy■^col 


C.3  Integrated  Coefficients 


The  integrated  coefficients,  col  to  co78.  shown  in  matrices  and  are  the 
results  of  the  product  of  functions  ^  and  rj  integrated  from  z  =  0  to  r  =  /. 


Functions  f  and  as  well  as  the  relations  between  con  (n  =  1 . 78)  and  functions 

T]  are  given  in  table  C-1. 


From  table  C-1.  we  arrive  at  the  following  values. 

kkl=0.9825022U5762 
klc2=l  .000’’773 11907 
kl =4. 730040744863 
k2=7. 853204624096 
col=0. 5 
co2=0 

co3=0. 2830976246421625 
co4=-  0.2828777404062893 
co5=-  79.9432918768421 
co6=-  79.88119926948301 
co7=0 . 3672702673507307 
co8=0 . 3608438510199932 
co9=-  5.754332755625593 
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Table  C-l:  Functions  ^  and  rj 
(Clamped-CIamped,  Nonnal  Mode  Functions) 
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colO=-  5.653644675810762 
col  1=0 

col2=0. 4244131815783876 
col3=0.5 

col4=0 . 353562193442563 
col5=0 . 3538370215454286 
col6=99 . 84162058130237 
col7=99. 91922866182124 
col8=-  0.2716466305565636 
c 019=0.2764844969574014 
c  o20=4 . 256116661561058 
co21=4 . 331915591047025 
co22=-  0.2122065907891938 
co23=0 

c024=0 .4999505291404635 
co25=0 . 06366823762540648 
co26=81 . 87642671677981 
co27=81 .94007023819137 
co28=0. 002279677482065308 
co29=0 . 3201666018578289 
co30=3 . 184351735522682 
co3l=3. 277293610756727 
co32=-  0.06064472050596255 
co33=0. 1517146400099545 
co34=0 . 5000494708595365 
co35=81 .94007023816478 
co36=82 . 00376323044337 
co37=-  0.3201552871859263 
co38=0 . 003620663312674233 
co39=5. 239125787449306 
co40=5 . 296229223177236 
co41=-  0.06069186036931498 
co42=-  0.1515968020106723 
co43=105455 . 1775540691 
co44=105455 .6457199359 
co45=-  121.5146304049034 
co46=-  70.38991464795974 
co47=5792 . 223435488106 
co48=5795 . 660646120006 
c 049=56 .51445896834048 
co50=-  141.3819822435741 
co51=105456 . 1775540691 
co52=-  121.4583153781003 
co53=-  70.29649790838054 
co54=5792 . 336747342393 


co55=5795 . 853^28-156008 
co56=56 . 5583883302294-1 
co57=-  141.2721697038976 
co58=0 . 5018655334042408 
co59=0 . 105674401288916 
co60=-  6.093554325721313 
co61=-  5.986930619661868 
co62=0. 2512142798792595 
c 063=0.3849690308126355 
co64=0. 4981344665957592 
co65='  5.986930619658262 
co66=-  5.88217259232853 
co67=-  0.2556882581559216 
co68=0 . 3782329253165705 
co69=338 . 7394246092603 
co70=339. 1865750592409 
co7l=5 . 519283205241623 
co72=-  8.457931242296016 
c  073=339 . 7394246092603 
co74=5 .617578386448839 
co75=-  8.309936176289431 
co76=0 . 5 
co77=0 
co78=0.5 


APPEP^DIX  D 


MATRICES  ME  AND  MF 


D.l  Matrix  Me 


M,  (ij)  =  M,  (;.i)  for  i=j  {tj=\ . 8) 

me(l  ,1)  =  ( ( (2''Co70-2“co65-2‘‘Co61+2*co59)*ge*kl*kkl**2+(-2'!'Co73-2* 

1  co69+-«>‘co66-2>'co6‘4+2*co58)‘’g€“kl*kkl  +  (2*co70-2*co65-*’2*co61-2*co 

2  59)  *-ge>"kl)  ‘'lkx‘2+  (2>‘Co70-2''Co65+2~co61*2*co59)*kl**3*kkl**2+(-2 

3  *co73-2“co69+2‘‘co64-4“co60-2*co58)*kl**3*kkl+(2*co70+2*co65+2“c 

4  061  +  2*0  059) *kl  '‘*3)*ll*xz+( ( (co73'2*co66+co64)*ge*kl**2»kkl**2+( 

5  -2*co70+2»co65-2*co61+2*co59)*ge*kl**2>«kkl+(co69+2+co60+co58)*g 

6  e*kl''*2)  *lk**2+(co69+2*co60+co58)*kl><*4»kk^**2+(-2*co70-2*co65- 

7  2*co61 -2*0059) *kl**4*kkl+(oo73+2*co66+co64)*kl**4) *xy+( ((co69-2 

8  *co60+oo58)'ge*kkl**2+(-2*oo70+2*co65+2»co61-2*co59)*ge*kkl+(oo 

9  73  - 2*oo66+c o64 ) *ge ) *lk**2+ (oo73-2*c 066+0 o64)»kl**2*kkl**2+( -2*0 
o70+2*oo65-2*oo6l+2*oo59)*kl**2*kkl+(co69+2*co60+co58)*kl**2)*l 

;  1**2+ ( (oo69-2*oo60+oo58)*kkl**2+(-2*co70+2*oo65+2*co61-2»oo59)* 

<  kkl+oo73-2*oo66+oo64)*la*lk**2 


meCl ,2)  =  (( (( (oo54-oo52-oo39+oo37)*ge*k2+(co48-co46-co31+co29)»ge 

1  *kl ) *kkl + ( -oo55+co53+oo40-oo38) *ge*k2+ ( -co47-oo45+co30+co28) *ge 

2  *kl)*kk2+((-co47  .o-.5-oo30+oo28)*ge*k2+(-oo55+oo53+oo40-co38)*g 

3  e*kl)  *kkl+(oo48-o>.  ♦co31-oo29)*ge*k2+(co54+co52-co39-co37)*ge* 

4  kl) *lk**2+( ( (oo48-oo46+oo3l-co29)*kl*k2**2+(co54+co52-co39-co37 

5  )*kl**2*k2)*kkl+(-oo47-co45-co30-co28)*kl*k2**2+(-co55-co53+co4 

6  0+CO38) *kl **2*k2) *kk2+ ( ( -co55+co53-co40+co38) *kl*k2**2+ (-co47-c 

7  o45-oo30-co28)*kl**2*k2)*kkl+(co54+co52+co39+co37)*kl*k2**2+(co 

8  48+co46+oo31+oo29) *kl**2»k2)*ll*xr+ ( ( ( (co55-co53-co40+co38) *g*+ 

9  kl*k2*kkl+(-co54-co52+oo39+oo37)*ge*kl»k2)*kk2+(-co48+co46-co31 
+co29)*ge*kl*k2*kkl+(co47+co45+co30+co28)*g«*kl*k2)*lk**2+((co4 

;  7+co45+co30+co28)*kl**2*k2**2*kkl+(-co48-co46-co31-co29)*kl’''*2* 

<  k2**2)*kk2+(-co54-co52-co39-co37)*kl**2*k2**2*kkl+(co55+co53+co 

=  40+CO38) *kl **2*k2**2) *xy+ ( ( ( (oo47-co45-co30+co28) *ge’fkkl+ ( -co48 
>  +co46+oo31-oo29) *g€) *kk2+(-oo54+co52+co39-co37)*ge*kkl+(co55-co 

?  53-oo40+oo38)*ge)*lk**2+((oo55-oo53-co40+co38)*kl*k2*kkl+(-co54 
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-co52+co39+co37)  »kl='k2)  *kk2+(-co48+co46-co31+co29)  “kl*k2>'kkl+  (c 

1  o47+co-i5+co30+co28)  “kl^k2)  •*ll*’^2+((  (co47-co45-co30+co28)‘‘kkl-co 

2  48+co46+co3l-co29)  »kk2+(-co54-'-co52+co39-co37)  >“kkl-*-co55-co53-co4 

3  0+co38)  xla*lk>‘*2 

raed  .3)  =  ( (((-co69+2»co60-co58)»‘ge*kkl**2+(2*co70-2*co65-2*co61+2 

1  *co59)  ■“gexkkl+('Co73+2>*'Co66-co64)*ge)’^lk*'*2+(-co73+2*co66-co64) 

2  *kl*x2xkklxx2+(2xco70-2xco65+2xco61-2*co59)*kl**2*kkl*(-co69-2* 

3  c  260-CO58) xklxx2)*ll*x2*( (co73-2*co66+co64)*ge*kl**2*kkl**2+(-2 

4  'i>co70+2»co65-2xco6l+2*co59)*ge*kl**2*kkl+(co69-*-2*co60+co58)*ge* 

5  klx*2)xlk*x2+(co69+2*co60+co58)*kl**4*kkl**2+(-2*co70*2*co65-2* 

6  co61-2xco59)xkl**4*kkl+(co73+2»co66+co64)*kl**4)*X2+(({-co70+co 

7  65+co61-co59)*gexkl*kkl**2+(co73+co69-2’'‘co66+co64-co58)*gexkl*k 

8  kl+(-co70+co65-co61+co59)*ge*kl)*lk»»2+(-co70+co65-co61+co59)*k 

9  lxx3xkklxx2+(co73+co69-co64+2xco60+co58)»klx*3*kkl+(-co70-co65- 
co61-co59)xklxx3)*ii»xy+(((co70-co65-co61+co59)*ge*kl*kkl*x2+(- 

;  co73-co69+2xco66-co64+co58)xge*kl*kkl+(co70-co65+co61-co59)*gex 

<  kl) xlkxx2+(co70-co65+co61-co59)*kl*x3*kkl**2+(-co73-co69+co64-2 
=  xco60-co58)  xklx»3xkkl+  (co70+co65-t'co61+co59)xkl*x3)  *11 


me(l ,4)  =  ( ( ( ( (-co47+co45+co30-co28)»ge>«kkl+(co48-co46-co31+co29) * 

1  ge)xkk2+(co54-co52-co39+co37)»ge*kkl+(-co55-<-co53+co40-co38)>*‘ge) 

2  *lkxx2-t'((-co55+co53+co40-co38)xkl*k2*kkl+(co54+co52-co39-co37)* 

3  klxk2)xkk2+(co48-co46-*'Co31-co29)»kl*k2»kkl+(-co47-co45-co30-co2 

4  8)*klxk2)xllx*2+(((co55-co53-co40+co38)’*‘ge’»klxk2xkkl  +  (-co54-co5 

5  2*co39+co37)»gexklxk2)*kk2+(-co48+co46-co31'*’Co29)*ge*kl»k2*kkl  + 

6  (co47+co45+co30+co28)“ge*kl*k2)*lk*x2+( (co47+co45+co30+co28) ^kl 

7  *x2><k2x*2*kkl+(-co48-co46-co31-co29)*kl*x2»k2*»2)*kk2+(-co54-co 

8  52-co39-co37) Xklx*2*k2xx2xkkl+(co55+co53+co40+co38)*kl*x2*k2x‘2 

9  )xxz+( ( ( (-co48+'co46+co31-co29)»ge*kl*kkl+(co47+co45-co30-co28) * 
ge*kl) xkk2+(co55-co53-co40+co38)*ge*kl»kkl+(-co54-co52+co39+co3 

;  7)xgexkl)xlkx*2+( (-co54-co52+co39+co37)*kl*x2*k2*kkl+(co55+co53 
<  -co40-co38)xkl**2xk2)xkk2+(co47+co45+co30+co28)*kl*x2*k2*kkl+(- 

=  co48-co46-co31-co29)*klx»2*k2)»ll*xy+((((co54-co52-co39+co37)xg 

>  exk2xkkl+(*co55+co53+co40-co38)*ge*k2)*kk2+(-co47+co45-co30+co2 

?  8)xge*k2*kkl  +  (co48-co46+co31-co29)*ge*k2)*lk>'‘*2+((co48-co46-*-co3 

l-co29)’i‘kl*k2**2*kkl+(-co47-co45-co30-co28)*kl*k2**2)*kk2+{-co5 

1  5+co53-co40-t-co38)*kl*k2**2*kkl  +  (co54+co52+co39i-co37)*kl*k2**2)* 

2  11 

me  (1.5)  =  (( (co68-co75)*kl*kkl+(co72+co63)*kl)*lk*ll*p2+((co62-co7 

1  l)*ge*kkl+(co74-co67)*ge)*lk**3*ll)*x2+(((-co72-co63)*kl**2*kkl 

2  +(co75+co68)*kl**2)*lk*p2+((co67-co74)’»ge*kl*kkl+(co71+co62)*ge 

3  *kl ) *lk**3) *xy+ ( (co72*co63) *kkl-co75+co68) *la*lk*p2 
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me  (1.6)  =  (((3«col9-3>‘Co21)-^kl*kkl+(3''co20+3*col8)*kl)*lk*ll*p2t-(( 

1  co7-co9)xge»kkl+(colO-co8)*ge)*lk**3*ll)*xz+(((-3*co20-3*col8)* 

2  kl'*''2^kkl-*-(3*co2l+3*col9)»kl**2)*lk*p2+(  (co8-col0)  “‘ge*kl*kkl+(c 

3  o9+co7)>«ge>-kl)xlk*>‘3)“xy+((3«co20*3*col8)*kkl-3>^co21+3*col9)»la 
*i  *lk>‘p2 

me(l  ,7)  =  ( ( (co68-co75)*ge>“kl«kkl+(co72+co63)*ge*kl)*lk*p2+((co62- 

1  co71)*ge*kkl+(co74-co67)»ge)*lk#ll**2+((co71+co62)*kl**2*kkl+(- 

2  co74-co67)»'kl**2)‘*'lk)»X2+((co67-co74)*ge*ki*kkl  +  (co71+co62)*ge* 

3  kl)‘'lk*ll*xy+((co63-co72)»ge*kkl+(co75*co68)*ge)*lk*ll*p2+((co7 

4  4-C067)  »kl*'kkl+(-co71-co62)*kl)*lk*ll 

me(l ,8)  =  (((3*col9-3*co21)*ge*kl*kkl+(3*co20+3*col8)*ge*kl)»lk*p2 

1  +( (co7-co9)  *ge^kkl+(col0-co8)*ge)»lk*ll*'i'2+((co9+co7)*kl<“‘2*kkl 

2  +(-co8-colO)  *kl*>‘2)»‘lk)''X2+((co8-colO)»ge*kl*kkl+(co9+co7)*ge*k 

3  1)  “Ik^ll^xy+C  (3»col8-3*co20)’<ge*kkl+(3>*co21-3*col9)*ge)*lk‘»ll*p 

4  2+(  (col0-co8)*kl‘^kkl  +  (-co9-co7)*kl)>!‘lk’*ll 

me  (2. 2)  =  (((2xco44-2»co35-2»-co2742>i<co25)>^ge*k2*kk2>‘*2+(-2-co51-2* 

1  co43+4''co36-2>‘co34+2*co24)*ge*k2»kk2+(2*co44-2>»co35+2*co27-2*co 

2  25)xge-k2)*lk-*24(2xco44-2*co35+2*co27-2*co25)*k2*»3*kk2*‘2+(-2 

3  >‘co51-2xco43+2“co34-4'“co26-2'«co24)*k2*<‘3*kk2+(2*co44+2*co3542’''c 

4  o27+2>‘Co25)»k2x*3)>‘ll‘‘Xz4(((coSl-2*co36+co34)*ge*k2**2*kk2*»2+( 

5  -2‘‘co44+2»co35-2‘'co27+2»co25)*ge’‘k2>«>‘2*kk2+(co43+2*co26+co24)’*‘g 

6  e*k2*>'2)»lk‘*24(co43+2‘‘Co26+co24)»k2**4*kk2‘“*2+(-2*co44-2*co35- 

7  2-co27-2«co25)’'k2*>'4‘>kk2+(co5l42»co36+co34)*k2*M)*xy+(((co43-2 

8  >'Co264co24)*ge>‘kk2>«»^2+(-2*co4442*co35+2*co27-2*co25)*ge*kk24(co 

9  51-2»co36+co34)*ge)*lk‘'»2+(co5l-2‘‘Co36+co34)*k2**2*kk2**24(-24c 
o44+2>‘Co35-2>»co27+2»co25)»k2>‘*‘2*kk2+(co43+2*co26+co24)*k2«*2)*l 

;  lx'‘24( (co43-2xco26+co24)*kk2“*2+(-2*co44+2*co35+2*co27-2*co25) * 

<  kk2+co51-2xco364co34)xla*lk*x2 

me (2 ,3)  =  ( ( ( ( (-co47+co454co30-co28)*ge*kkl+(co48-co46-co31+co29)“ 

1  ge) xkk2+(co54-co52-co39+co37)*ge*kkl+(-co55+co53+co40-co38) *ge) 

2  *lkx*24((-co55+co53+co40-co38)*kl*k2*kkl+(co54+co52-co39-co37)* 

3  kl*k2)*kk2+(co48-co46+co31-co29)*kl*k2*kkl+(-co47-co45-co30-co2 

4  8)»kl*k2)*ll**2+(((co55-co53-co40+co38)*ge*kl*k2*kkl+(-co54-co5 

5  24co394co37)*g«*kl*k2)*kk2+(-co484co46-co3l4co29)*g**kl*k2*kkl+ 

6  (co474co45+co304co28)*ge*kl’'‘k2)*lk**2+(  (co47+co45+co30+co28)*kl 

7  **2*k2**2*kkl4(-co48-co46-co31-co29)*kl**2*k2**2)’'‘kk24{-co54-co 

8  52-co39-co37)*kl»*2*k2**2*kkl+(co55+co53+co40+co38)*kl**2*k2**2 

9  )*xz+((((-co54+co62+co39-co37)*g«*k2*kkl+(co66-co53-co404co38)* 
ge*k2)*kk2+(co47-co45+co30-co28)*g**k2*kkl+(-co484co46-co31+co2 

:  9)*ge*k2)*lk**2+((-co48+co46-co3l4co29)*kl*k2**2*kkl+{co47+co45 

<  4co304co28)*kl*k2**2)>*kk2+(co55-co534co40-co38)‘''kl»k2**2*kkl  +  (- 
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=  co54-co52-co39-co37)*kl’^k2**''2)»‘ll*X)'+  ( ( ( (co48-co46-co31+co29)*g 

>  e*kl*kkl  +  (-co47-co45+co30-»co28)*ge*kl)*‘kk2+(-co55+co53+co40-co3 

?  8)«ge*kl'*kkl+(co54+co52-co39-co37)>»ge*kl)*lk*»2+(  (co54+co52-co3 

9-C037)  >‘kl>“>‘2*k2“^kkl-»-(-co55-co53+co40+co38)*kl**2*k2)*kk2+  (-co4 

1  7-co45-co30-co28)  >*klx*2*k2*kkl+(co48+co46+co31+co29)*kl**2*k2)* 

2  11 

me(2,4)  =  ( ( ( (-co43+2*co26*co24)*ge*kk2*»2+(2*co44-2*co35-2*co27+2 

1  *co25)-ge*kk2+(-co51+2»co36*co34)»ge)*lk»*2+(-co51+2*co36*co34) 

2  *k2**2‘»kk2**2+(2*co44-2*co35+2*co27-2*co25)*k2**2*kk2+(-co43-2» 

3  co26-co24)*k2**2)*ll**2+((co51-2*co36+co34)*ge*k2**2*kk2**2+(-2 

4  *co44+2‘^co35-2*co27'f2*co2S)*ge»k2**2*kk2+(co43+2*co26+co24)*ge* 

5  k2'i""2)*lk**2+(co43+2*co26+co24)*k2*»4*kk2**2+(-2*co44-2»co35-2* 

6  co27-2*co25)*k2‘<*4>'kk2+(co51+2*co36+co34)*k2**4)*x2+((  (-C044+CO 

7  35+co27-co25)»ge*k2*kk2-^*2+(co51+co43-2*co36+co34-co24)*ge*k2*k 

8  k2+{-co44+co35-co27+co25)’i-ge*k2)*lk»-!‘2+(-co44+co35-co27+co25)*k 

9  2»'»3»kk2x>i'2+(co51-*'Co43-co34+2*co26+co24)*k2**3*kk2+(-co44-co35- 
C027-C025)  *k2»>'3)  *ll^xy+  ( ( (co44-co35-co27+co25)  »ge*k2*kk2»>“2+  ( - 

:  co51-co43+2"Co36-co34+co24)*ge’*'k2»kk2+(co44*co35-*-co27-co25)»‘ge'‘ 

<  k2)*lk**2+(co44-co35+co27-co2S)*k2*>*3*kk2*>»2+(-co51-co43+co34-2 

=  *co26-co24)*k2**3*‘kk2+(co44+co35+co27+co25)*k2**3)*ll 


me  (2 ,5)  =  ( ( (co42-co57) *k2*kk2+(co50+co33) *k2) »lk*ll*p2+ ( (co32-co4 

1  9)*ge»kk2+(co56-co41)*ge)*lk*«3»ll)><xz+(((-co50-co33)*k2*»2»kk2 

2  +(co57*co42)*k2»*2)*lk*p2+(  (co41-co56)‘i'ge»k2*kk2+(co49+co32)»ge 

3  *k2)  *lk>'*3)  >«xy+  ( (co50-co33)  *kk2-co57+co42)  *la>i‘lk’*p2 

me  (2. 6)  =  (((3»col5-3‘‘Col7)*k2*kk2+(3*col6+3*col4)*k2)*lk‘fll‘»p2+(( 

1  co3-co5)»gexkk2+(co6-co4)'^ge)‘i‘lk*‘*3’'ll)*xz+(((-3*col6-3*col4)*k 

2  2»*2-i<kk2+(3'^col7+3'‘Col5)*k2»*2)*lk*p2+((co4-co6)*ge*k2*kk2+(':o5 

3  -•■C03)  ’“ge»k2)  *lk’i‘’^3)  *xy+(  (3*col6-3*col4)  *kk2-3*col7+3*col5)>^la*l 

4  k*p2 

me(2,7)  =  ( ( (co42-co57)*ge*k2*kk2+(co50+co33)*ge*k2)*lk*p2+( (co32- 

1  co49)*ge*kk2+(co56-co41)*ge)*lk*ll**2+((co49+co32)*k2**2*kk2+(- 

2  co56-co41)*k2**2)*lk)*xz+((co41-co56)*ge*k2*kk2+(co49+co32)*ge* 

3  k2)*lk*ll*xy+((co33-co50)*ge*kk2+(co57-co42)*ge)»lk*ll*p2+((co5 

4  6-co41)*k2*kk2+(-co49-co32)*k2)*lk*ll 

me (2, 8)  =  ( ( (3*col5-3*col7)*ge*k2*kk2+(3*col6+3*col4)*ge*k2)*lk*p2 

1  +  ( (co3-co5)  *ge*kk2+  (co6-co4)  *ge)  ■'‘lk*ll**2+  ( (co5+co3)  ’»k2**2*kk2+ 

2  ( “Co6-co4) *k2**2) *lk) *X2+ ( (co4-co6) *g«*k2*kk2+ (co6+co3) *ge*k2) * 

3  lk*ll*xy+((3*col4-3*col6)*g**kk2+(3*cel7-3*col5)*ge)*lk»ll’''p2+( 

4  (co6-co4)*k2-*'kk2+(-co5-co3)*k2)*lk*ll 
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me (3 ,3)  =  ( { (-2‘'Co70+2>'Co65-»2*co61-2*co59)*‘ge’'kl*kkl**2-*'(2“-co73+2* 

1  co69-4*co66+2>‘co6‘;-2“co58)*ge*kl*'kkl+  (-2’*co70+2*co65-2*co61+2*c 

2  o59)  “ge^kl)>'lk*'<2+(-2‘‘co70+2*co65-2*co61+2*co59)*kl*»3*kkl**2+  ( 

3  2*co73+2‘^c 069-2*0 o64+4“co60+2*co58)*kl**3*kkl+(-2*co70-2*co65 *2 

4  *co61-2*co59)*kl**3)*ll*X2+(((co69-2*co60+co58)*ge*kkl**2+(-2*c 

5  o70+2*co65+2>^co61-2^co59)»ge»kkl+(co73-2*co66+co64)*ge)*lk**2+( 

6  co73-2*co66+co64) *kl**2*kkl*»2+(-2*co70+2*co65-2*co61+2*co59)*k 

7  1  **2>‘kkl+  (co69+2>^co60+co58)  *kl**2)  *ll**2*xy+  ( (co73-2*co66+co64) 

8  *ge*kl**2*kkl**2+(-2*co70-''2«co65-2*co61+2*co59)*ge*kl**2*kkl  +  (c 

9  o69+2*co60+co58)=»ge*kl**2)*lk**2+(co69+2*co6O+co58)*kl**4*kkl** 
2+{-2*co70-2*co65-2“»co61-2“^co59)»kl*‘'‘4*kkl  +  (co73+2’'‘co66+co64)*k 

:  1**4 


me(3.4)  =  ((( ((-co54+co52+co39-co37)*ge*k2+{-co48+co46+co31-co29)* 

1  ge*kl ) *kkl + (co55-<:o53-co40+co38) *ge*k2+ (co47+co45-co30-co28) *ge 

2  *kl ) *kk2+ ( (co47-co45+co30-co28) *ge*k2+(co55-co53-co40+co38) *ge* 

3  kl) *kkl+(-co48+co46-co31+co29)*ge*k2+(-co54-co52+co39+co37)*ge* 

4  kl)*lk**2+( ( (-co48+co46-co31+co29)*kl*k2**2+(-co54-co52+co39*co 

5  37) *kl>'*2*k2) *kkl+(co47+co45+co30+co28)*kl*k2**2+(co55+co53-co4 

6  0-CO38)  *kl**2>‘k2)  *kk2+(  (co55-co53+co40-co38)*kl*k2**2-*-(co47+co4 

7  5+co30+co28)>‘kl **2*k2)*kkl+{-co54-co52-co39-co37)*kl*k2**2+(-co 

8  48-co-,6-co31-co29)*kl**2>'k2)*ll*xz+((((co47-co45-co30+co28)*ge* 

9  kkl+(-co48fco46+co31-co29)*ge)*kk2+(-co54+co52+co39-co37)*ge*kk 
l  +  (co55-co53-co40+co38)*ge)*lk**2+((coB5-co53-co40+co38)’i‘kl*k2* 

:  kkl*(-co54-co52*co39*co37)*kl*k2)*kk2+(-co48+co46-co31+co29)*kl 

<  *k2‘kkl+(co47+co45+co30*co28)“kl*k2)*ll**2*xy+( ( (co55-co53-co40 

=  +0038) *ge*kl *k2*kkl * ( -co54-co52+co39+co37) *ge»kl*k2) *kk2+ ( -co48 

>  +co46-co31*co29)*ge*kl»k2*kkl+(co47+co45+co30+co28)*ge»kl*k2)*l 

?  k**2+( (co47->co45+co30+co28)*kl**2*k2**2*kkl+(-co48-co46-co31-co 

29)*kl*'‘2*k2**2)>‘kk2+(-co54-co52-co39-co37)*kl**2*k2**2*kkl  +  (co 
1  55+co53+co40+co38) *kl**2*k2**2 

me (3 ,5)  =  ( ( (-co72-co63) *kl**2*kkl+(co75+co68)*kl**2) *lk*p2+( (co67 

1  -co74)*ge*kl*kkl+(co71+co62)*g€*kl)*lk**3)*xz+( ((co75-co68)*kl* 

2  kkl+(-co72-co63)*kl)*lk*ll*p2+((co71-co62)*ge*kkl+(co67-co74)*g 

3  e)*lk**3*ll)*xy 

me (3, 6)  =  (( (-3*co20-3*col8)*kl**2*kkl+(3*co21+3*col9)*kl**2)*lk*p 

1  2+( (co8-col0)*ge*kl»kkl+(co9+co7)*ge*kl)*lk**3)*xz+(({3*co21-3* 

2  col9)*kl*kkl+(-3»co20-3»col8)*kl)*lk*ll*p2+((co9-co7)*ge*kkl+(c 

3  o8-col0)*ge)*lk**3*ll)*xy 

me  (3, 7)  =  (( (co72-co63)*ge*kkl+(co68-co75)*g*)*lk*ll*p2+(((co67-co 

1  74)*ge-co74+co67)*kl*kkl+((co71+co62)*g«+co7l+co62)*kl)*lk*ll)* 

2  x2+( (co71 -co62)*ge*kkl-^(co67-co74)*ge)*lk*ll**2*xy+{(co68-co75) 


92 


3  ~gexkl'‘kkl+(co72+co63)''ge"kl)>-lk*-p2+((co7l+co62)*klxx2»kkl  +  (-co 

4  74-co67)-kl>^^2)-lk 

me  (3, 8)  =  (( (3*co20-3>‘Col8)*ge»kkl  +  (3’^col9-3*co21)*ge)*lk>‘ll>"p2+(  ( 

1  (co8-colO)~ge-t-co8-colO)’^kl*kkl+((co9+co7)*ge+co9+co7)>i‘kl)*lk*ll 

2  )>»xz+(  (co9-co7)»ge*kkl+(co8-colO)*ge)*lk*ll**2'i‘xy+(  (3*col9*3*co 

3  21)>'ge*kl*kkl+(3*co20+3*col8)*ge^kl)*lk*p2+( (co9+co7)*kl**2*kkl 

4  +(-co8-col0)*kl**2)*lk 

me  (4, 4)  =  (( (-2>‘co44+2*co35+2*co27-2*co25)*ge*k2*kk2**2+(2*co5l+2* 

1  co43-4-“co36+2*co34-2»co24)  *ge*k2*kk2-*-(-2*co44+2*co35-2>*co27+2>^c 

2  o25)*ge*k2)*lk*>‘2+(-2*co44+2^co35-2*co27+2*co25)*k2**3*kk2*»2*( 

3  2^co5l+2*co43-2'‘co34+4*co26+2*co24)*k2**3*kk2+(-2*co44-2»co35-2 

4  *co27-2*co25)‘<k2*>'3)‘'ll*xz+((  (co43-2*co26+co24)*ge*kk2**2+(-2*c 

5  044+2-0 o35+2*co27- 2* co25)  ‘‘ge*kk2+(co51-2*co36+co34)*ge)*lk+»2+( 

6  co5l-2*'co36+co34)>'k2»*2*kk2*»2+('2*co44+2*co35-2*co27+2*co25)  ‘■k 

7  2*“2*kk2+(co43+2«‘co26+co24)“k2»+2)+ll**2+xy+((co51-2*co36+co34) 

8  >‘ge-k2**2“^kk2-+2+(-2*co44+2«co35-2‘co27+2+co25)*ge*‘k2*«2*kk2+(c 

9  o43+2>'co26+co24)*ge+k2**2)*lk*»2+(co43+2*co26+co24)*k2**4*kk2*‘>‘ 
2+(-2>^co44-2»co35-2*co27-2*co25)>‘k2+*4*kk2+(co51+2+co36+co34)»k 

;  2-'‘4 


me  (4 ,5)  =  ( ( (-co50-co33)*k2**2*kk2+(co57+co42)*k2*’^2)'“lk*p2+(  (co41 

1  -co56)  ’"ge*k2*kk2+(co49+co32)>“ge*k2)‘<lk*+3)*X2+(  ((co57-co42)-k2- 

2  kk2+(-co50-co33)-k2)*lk‘'ll*p2+((co49-co32)*ge+kk2+(co41-co56)-g 

3  e)-lk->‘3'll)-xy 

me  (4 .6)  =  ( ( (-3»col6-3*col4)>‘k2**2*kk2+(3*col7+3*col5)*k2’'  +  2)*lk>-p 

1  2+((co4'Co6)»ge»k2*kk2+(co5+co3)*ge*k2)*lk**3)*xz+(  ( (3>'col7-3’'C 

2  0l5)*k2-kk2+(-3‘>col6-3-col4)*k2)*lk*ll»p2+((co5-co3)*ge+kk2+(co 

3  4-co6)  «ge) -lk-^3’^11)  “xy 

me (4 ,7)  =  (( (co50-co33)»ge*kk2+(co42-co57)+ge)*lk*ll»p2+(((co41-co 

1  56)  *ge-co56+co41 )  *k2=»kk2+  ( (co49+co32)  ‘ge+co49+co32)  *k2)  *lk*ll) 

2  xz+( (co49-co32)*ge*kk2+(co41-co56)*ge)*lk*ll**2*xy+((co42-co57) 

3  *ge-k2+kk2+(co50+co33)+ge»k2)*lk*p2+((co49+co32)*k2**2*kk2+{-co 

4  56-co41)>^k2*+2)*lk 

me (4, 8)  =  (( (3+col6-3*col4)*ge*kk2+(3*col5-3*col7)*ge)*lk*ll*p2+(( 

1  (co4-co6)*ge-co6+co4)*k2*kk2+((co5+co3)*g«+co5+co3)*k2)*lk*ll)* 

2  Xz+((co5-co3)*ge*kk2+(co4-co6)*g«)*lk*ll*’*2*xy+((3*col5-3*col7) 

3  *g€*k2*kk2+ (3*col6+3*col4) *ge*k2) *lk*p2+ ( (coB+co3) *k2**2*kk2+ ( - 

4  co6-co4)*k2**2)*lk 


me(5,5)  =  (co78*lk>“*2*p2**2+co76*g«»lk**4)*xy+co78*la*p2**2 
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me  (5 ,6)  =  (3*co23-lk-^2-p2'‘-2+coll-'ge»lk'">4)»xy+3^co23*la»p2*>‘2 

me(5.7)  =  (coTT-ge-co??) -lk'<''2'‘p2-'Xz*co76>‘ge>'lk“>*2“ll*xy 

me(5,8)  =  (3“co22’ge-col2)>‘lk‘“'2’'p2>‘X2+coll’!‘ge*lk**2>‘ll>*xy 

me(6,6)  =  (9‘‘Col3'lk>'>'2>‘p2"»2+col“ge*lk*>‘4)*xy+9*col3«la“i'p2**2 

me(6,7)  =  (col2>'ge-3«co22)*lk*“2*p2'‘xz+coll*ge*lk**2*ll‘'xy 

me(6.8)  =  (3>‘Co2“ge-3*co2)  »lk>»>'2*p2“XZ+col*ge*lk**2»ll*xy 

me (7 , 7)  =  -Imda*-  (-co76“xy-co76)*(co78»ge*p2**2-co76*lmda+co76*ge*l 
1  1"''2)  "X\’-'co78“ge‘'p2>'>-2-co76»lmda*co76*ge*ll**2+co76»lk’i'’^2 

me(7  .8)  =  -Imda'‘(-coll“xy-coll)  +  (3>'co23’>ge*p2“*2-coll*lmda+collxge 
1  *11"‘2)  >‘xy+3'-co23>‘ge>‘p2*>' 2-col  l~lmda+coll*ge«‘ll*»2+coll»lk*>‘2 

m€(8.8)  =  -Imda'- (-col“xy-col)  +  (9»col3-’ge“p2**2-col*ljnda+C'  »e»ll» 
1  '2)  '‘X^'+9-col3“ge>‘p2“*2-col»lmda+col*ge*ll**2+col*lk'‘*2 


D.2  Matrix  Mf 

M^(ij')  =  for  j=5  (ij=l . S) 

mf  (1 , 1)  =  -  ( (-co69*2-'co60-co58)>'kkl*“‘2+(2“co70-2“co65-2*co6l  +  2»co5 
1  9)’'kkl-co73+2>‘Co66-co6‘i)-la 

mf(l,2)  =  -  ( ( (-co47*co.;5*co30-co28)  >‘kkl+co48-co46-co3l+co29)  «kk2+( 
1  co54-co52-co39+co37)»kkl -co55+co53+co40-co38)“la 

mf(1.3)  =  0 

mf  (1 ,4)  =  0 

mfd.S)  =  0 

infCl.e)  =  0 

mf(1.7)  =  0 


mf(1.8)  =  0 
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mf  (2.2)  =  -( ( -c 0-43+2' c o26- c o2-<) kk2’“‘2+ (2»c 044 -2*co35-2“Co27+2>^co2 
1  5) 'k^ . -:o51+2>‘Co36-co34) 'la 

';.3)  =  0 

ml(2.4)  =  0 

mf(2.5)  =  0 

mf(2.6)  =  0 

mf(2.7)  =  0 

mf(2.8)  =  0 

ml  (3 .3)  =  -  ( (-co69+2>-co60-co58)+kkl"2+(2*co70-2+co65-2-co61+2'Co5 
1  9) 'kkl -co73+2«co66-co64) 'la 

mf(3.4)  =  -  ( ( (-co47+co45+co30-co28)’»kkl+co48-co46-co31+co29) 'kk2+ ( 
1  co54-co52-co39+co37)»kkl -co55+co53+co40-co38)*la 

ml(3.5)  =  0 

mf(3.6)  =  0 

ml(3.7)  =  0 

mf(3.8)  =  0 

mf  (4 .4)  =  -  ( (  -  co43+2«co26-co24)  ‘■kk2"2+(2'co44-2'co35-2'co27+2'co2 

1  5) -kk2-co51+2'co36-co34) 'la 

ml(4.5)  =  0 

ml(4,6)  =  0 

mf(4.7)  =  0 

ml(4,8)  =  0 

ml(5,5)  =  co76»la 

ml(5,6)  =  coll*la 


ml(5.7)  =  0 
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mf(5.8)  =  0 

mf(6,6)  =  col'la 

raf(6.7)  =  0 

mf(6.8)  =  0 

mf(7,7)  =  co76>'xy+co76 

mf(7.8)  =  coll>‘xy+col  1 

mf(8.8)  =  col-xy+col 


D.3  Integrated  Coefficients 


The  integrated  coefficients,  col  to  co78.  shown  in  matrices  > 
results  of  the  product  of  functions  ‘  and  rj  integrated  from  r  =  0  to  z  =  /. 


Functions  f  and  t]  as  well  as  the  relations  between  con  (n  =  1 
f.  rj  are  given  in  table  D-1. 


From  table  D-1.  we  arrive  at  the  following  values. 


kkl=0. 7340955137589 
kk2=l .018467318759 
kl=l. 875104068712 
k2=4. 694091 132974 
p2=0.5*pi 


col=0.5 

co2=0 . 1061032953945969 
co3=0. 4989995345474672 
co4=0 . 1 108751276626285 
co5=-  5.798956073142442 
co6=-  5.691469236183133 
co7=0. 03006018182478657 


j  and  My  are  the 


78)  and  functions 


sh(tJ>  sln<t1>  cosft^ 


Table  D-1:  Functions  ^  and  fj 
{Clamped -Free,  Normal  Mode  Functions) 
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co8=0 . 3-;78402254263256 
co9=-  0.2432880298097214 
col0=-  0.04890990761714533 
coll=0 

col2=0 . 3183098861837907 
col3=0.5 

col4=0 . 101726459217963 
col5=0 . 5009446645167207 
col6=-  5.926688707716224 
col7=-  5.821560750035544 
col8=-  0.3408655849552808 
CO 19=0. 07554528409550446 
co20=-  0.6562209160646115 
co2l=-  0.6114155741352192 
co22=-  0.3183098861837907 
co23=0 

co24=0 . 4980514051656076 
co25=0. 1064812333022709 
co26=-  5.714131203766796 
co27=-  5.606622101392816 
co28=0 . 03475931115663723 
co29=0 . 3486792690147293 
co30=-  0.2341827624244058 
co3l=-  0.03869781747047025 
co32=0 . 004389358274729287 
co33=0 . 3201612686560902 
co34=0 . 5019485948343924 
co35=-  6.032689670886772 
co36=-  5.927129326459122 
co37=-  0.3424954405873828 
co38=0 . 07770619523077635 
co39=-  0.6627446803520218 
co40=-  0.6173170841492832 
co41=-  0.3201345481538806 
co42=0 .00146882275175669 
co43=317 .6223970917853 
co44=318 . 0691431037289 
co45=10. 78126487405295 
co46=0 .6337238477409987 
co47=25. 64535716383292 
c 048=28 . 36465756278566 
co49=10. 47070965266181 
co50=3 . 31 1674398369181 
co51=318. 6223970917853 
c  o52= 10 . 85325055404595 


98 


co53=0. 8185862051268238 
co54=25. 73613685358632 
co55=28 .60744933975992 
co56=10 . 53306602928041 
co57=3 . 503841701282994 
co58=0 . 5762267152276445 
co59=0. 2427119121868963 
co60=l . 103441151283322 
co58=0. 5762267152276445 
co59=0 . 2427119121868963 
co60=l . 103441151283322 
co61=l .343334932063346 
co62=0 . 5357954242210862 
co63=0 . 3590249035362438 
co64=0. 4237732847723555 
co65=0 . 276727465644721 
co66=0 . 5946402690834619 
co67=0. 2080415754963285 
c  o68=0 . 4488420127305808 
co69=2 . 333604078798335 
co70=2 . 703413018359166 
co7l=l .045891914197238 
co72=0 . 5225168922652508 
co73=3 . 333604078798335 
co74*1 .260362987158822 
co75=0. 876155731545073 
co76=0.5 

co77=0 . 3183098861837907 
co78=0.5 
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